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SUMMARY 


A method has been developed which makes possible the de- 
termination of the transition time for an airplane to maneuver 
into steady turning flight from a level straight-line course and 
the control deflections necessary to execute this maneuver. Ac- 
celerations encountered are easily obtained from data necessary 
to compute the maneuvering time. The method developed is 
based on the six fundamental stability equations. The major 
assumptions are angles of bank, #, and no sideslip during turn, 
ie., v = 0. Only the use of conventional wind-tunnel data is 
necessary. 

Calculations are presented for a typical twin-engined attack 
bomber. The effect of mass distribution is shown by considering 
a modified single-engined airplane of the same weight and horse- 
power as that of the twin-engined design. The bank assumption 
chosen produced a helix angle of 0.090 and maximum bank of 
75°. Calculations were carried out for sea level and for 15,000 
ft. altitude. 

From the standpoint of practical application, the method de- 
veloped presents not only a means of estimating the maneuvering 
characteristics of a particular airplane but offers a basis for com- 
parison of similar airplanes, since control deflections necessary 
to effect a given maneuver are available as a function of time. 
This is illustrated in the calculations worked out for the two air- 
planes mentioned above. The time to maneuver from straight 
and level flight into a steady turn is substantially the same at sea 
level and 15,000 ft., but at high altitude the control deflections 
required for the twin-engined airplane exceed those necessary 
for the modified single-engined one. As suspected, this shows 
that an increase in spanwise mass distribution is detrimental to 
the maneuvering characteristics of an airplane 


INTRODUCTION 


P | ‘HE IMPROVEMENT of the flying qualities of modern 
aircraft is impeded by the absence of quantitative 
maneuverability and control criteria.'| Wing loading, 
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span loading, power loading, altitude, and speed are all 
factors influencing the behavior of the aircraft. To get 
a concise maneuverability criteria that takes into ac- 
count all these variables is extremely difficult, and no 
attempt has been made to do so in this report. 

As a special definition of maneuverability, the air- 
plane was arbitrarily assigned a flight path, and equa- 
tions were worked out to give the time necessary to 
execute the maneuver as well as the control motions. 

The flight path chosen was that executed by the air- 
plane in going from straight and level flight into a steady 
turn. For necessary simplification, application of con- 
stant power was maintained throughout the turn. 
This phase of maneuverability is of great importance 
in the case of a night fighter plane when it is in pursuit 
only by instrument and when a rapid turn at full power 
may be necessary. The basic assumption for this desig- 
nated flight path was the rate of change of angle of 
bank with time. This relationship was taken as an 
exponential function. 

Calculations are worked out for a typical, modern, 
twin-engined bomber, and the effect on maneuver- 
ability of engines mounted in the wing is shown by 
considering the airplane single-motored, the motor 
being mounted in the fuselage, having the same weight, 
and developing the same horsepower as the two engines 
installed in the wing of the initial airplane. 

Methods presented in this paper are based on the 
interpretation of wind-tunnel data that should always 
be available, at least in preliminary form, during the 
design stages of aircraft on which rational analysis is 
important. All curves have been calculated for a con- 
ventional airplane of about 20,000 Ibs. gross weight. 
Although the above magnitude is of no particular sig- 
nificance, the general trends and relationships should 
hold for airplanes of twice and half this size. 

Outside of general specifications and speed conditions, 
the only factors necessary are: 
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change in side force per unit mass due to 
unit rudder deflection 


In, = change in rolling moment per unit moment 
of inertia due to unit aileron deflection 

No, = change in yawing moment per unit moment 
of inertia due to unit rudder deflection 

No, = change in yawing moment per unit moment 
of inertia due to unit aileron deflection 

LI, = change in rolling moment due to yawing 
velocity 

L, = change in rolling moment due to rolling 
velocity 

N, = change in yawing moment due to yawing 
velocity 

N, = change in yawing moment due to rolling 
velocity 


Basic EQUATIONS AND ASSUMPTIONS 


Because of the bilateral symmetry of the airplane it 
is customary to divide the motions into two independ- 
ent groups—the lateral and the longitudinal—each con- 
sisting of three degrees of freedom: (1) lateral motions 
—(a) rolling, (b) yawing, (c) sideslipping; and (2) 
longitudinal motions—(a) pitching, (b) vertical trans- 
lation, (c) forward translation. 

Presumably, the reactions to small increments of 
longitudinal speed or displacement do not sensibly 
influence the lateral motions and the two groups may 
be treated independently. 


The Lateral Motions 

If the flight path is assumed to be horizontal (or 
nearly so) and the main forward velocity U» to be sub- 
stantially constant, the equations of motion in a lateral 
disturbance may be written :* 
= gsin® —rU,+vY,+ Vo 
vL, + pl, + rL, + Lo 
vN, + pN, + rN, + No 


(In sideslipping) dv/dt 
(In rolling) dp/dt = 
(In yawing) dr/dt = 

The quantities that arise in the consideration of the 
lateral motions are defined in N.A.C.A. T.R. No. 560 
and are in keeping with standard N.A.C.A. conven- 
tion. 

It has been found convenient to transform all sta- 
bility derivatives and disturbing effects into terms of 
accelerations of the airplane rather than retaining them 
as moments and forces. This transformation is ac- 
complished by dividing out the appropriate moments 
of inertia and the mass of the airplane. For ex- 
ample, 


(OL/0p)/mK,? may be written simply as L,; similarly, 
(ON/dv)/mK,? = N, and (OY/dr)/m = Y,. 


On the basis of these definitions the above equations 
for sideslipping, rolling, and yawing were written. 
Since the axes chosen will ordinarily lie near the axes 
of the principal moments of inertia of the airplane, 
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terms involving the products of inertia have been 


neglected. 


Solution of the Lateral Equations 


From the three lateral equations, the unknowns are 


, r, v, Yo, Io; No. 
One of the unknowns can be eliminated by letting 


Yo = 5, Yo, 

Io = baLog 

No = 5,No, + 5aNoa 

For the two other unknowns: (1) specify ®; (2) 
specify v. 

Taking v = O (as suggested in N.A.C.A. T.R. No. 


560) and assuming 


= Lowe — p-~(nt+m)t 
ean(1o ee Sa ) 





n n+m 


the three equations of lateral motion may now be 


written : 
rU) — br¥Yo, = gsin ® (1) 
rL, + baLog bd (dp/dt) ime pLp (2) 
(dr/dt) — rN, — 5,No, — 5gNoq = PNp (3) 
where 
Yo, = acceleration sideways due to unit rudder de- 
flection 
In, = rolling moment due to unit aileron deflection 
No, = yawing moment due to unit rudder deflection 
Nog = yawing moment due to unit aileron deflection 


The three unknowns in Eqs. (1), (2), and (8) are: r, 
6,, 5g. Multiplying Eq. (1) by No,/Yo,, and Eq. (2) by 
Noq/ La» 


6, = (—g/ Yo,) sin ® + (r/ Yo,) Uy (4) 
1 dp 1 1 ' 
ws. geen ae 5 
Ly, dt p a ” aa (5) 


Substituting 6, and 6, in Eq. (3) and recalling that ® = 


e(? —-e* 1i- stn, 





n n+m 





= + Ar = Ce~™ + Coe" t™™ — g mn ® (6) 
where 
C, = N,k a Ne kn — Nog L,k 
Log 0a 
C, = — N,k+ Nog k(n + m) + a L,k 
0g 


da 


so’ L, - No uy, — N, 
Lo Y, 


(6) 


Solving Eq. 
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aes e7 +m)t__ 


The boundary conditions are r = 0 at ¢ = O, there- 


fore, 


CQ is"? 


od ae — @~*) 4 
Y rt ) 


r 


Cy (e~ “+ m)t 
A —(n + m) 


—e~“*) — 


t 
greena f e“' sin @dt (7) 
Yo, 0 


The integral term in this equation must be evalu- 
ated: 

Let the true integral= fo‘ e4* sin ®(x) dx 

Let the approximate integral = sin ®(f) fo‘ e4* dx 


R = true/approximate 
Integrating the approximate integral 


_A S' 4-9 sin B(x) dx 
sin &(t) (1 — e~“*) 





For t small 


}3 
sin 6() = 6(¢) + | oe 
®(t) less than */2 
t less than '/, 


For simplicity let, 
LL=A fo e4*~- sin &(x) dx = A fo! e4* —? &(x)dx 


Integrating by parts, recalling again that 
ee at 


P(3 : 
(x) n+m 
—(n+m}t .~Af 
L, st pene. e ) 
A — (n+ m) 
e7(™ +m)t _ e7 At 
Peccesy, 


1 — e~4! (8) 





Therefore, for values of ¢ up to approximately '/, 
sec., 


Jo! e*! sin &(t) dt = Rsin &(t) fo e“* dt 
where R is given in Eq. (8). 
For t large 


The integral can be evaluated by an asymptotic 
evaluation. 
I, = e4* fo! e** sin (x) dx 
Expressing sin ®(x) in a Taylor’s series and sub- 
stituting in the above integral 


sin &(¢) + cos P(t) ie @-% 


sin ®(x) = 7 ~ 
t ! 


t A At] 
e** sin &(x) dx = sin ®(t) ar Sos + 
0 4 


At _ 
cos ®(f) zl3 _ a] 
dt LA A? 


Therefore 


sin P(t) 


=e {ih = oh) 


E cy (0) (*) (: — ¢~ 4! a) 
A /\dt 1—e* 


Summarizing 


2 = 





For ¢ less than '/¢ sec., 
. : . sin &(t 
oe f e“' sin b(t) dt = R oo (1— 
0 


For ¢ greater than '/2 sec., 


t 
g<* f e“* sin &(t) dt = (1 — e~4*) x 
0 


is - (ota) (ere guy 
A dt 1—e~* 


Therefore, by solution of the lateral equations of 
motion through assuming a bank ®(¢) in a horizontal 
turn in which there is no sideslip, the rudder deflection, 
5,, the aileron deflection, 5,, and the angular velocity 
of yaw, 7, can be found in terms of the airplane accelera- 
tions for unit control movement, Yo,, Lo,, No, No, and 
damping derivatives L,, L,, N,, N,. 


ee 


sin ®(t) 





The Longitudinal Motions 


The equations of motion in a longitudinal disturbance 
may be written, 


du/dt = uX, + wX,, — g0+ Xo 
dw/dt = uZ, + wZ, + qUo + Zo 
dq/dt = uM, + wM,, + gM, + Mo 


As before, the quantities that arise in the considera- 
tion of the longitudinal motions are defined in N.A.C.A. 
T.R. No. 560. 


Solution of the Longitudinal Equations 


From the three longitudinal equations, the unknowns 
are: u, w, 0, Xo, Zo, Mo. 

As in the case concerning the lateral equations, there 
are three equations and six unknowns. One of the un- 
knowns can be eliminated by letting Mp = 6,Mo, and 
Zo = 5 Zo, 

For the other two unknowns: (1) Xo = 0, i.e., no 
change in thrust; (2) @ = w/U); and (3) steady state, 
i.e., du/dt = 0, dw/dt = 0.* 


* Only the steady state condition of the turn will be considered, 
for that will determine maximum elevator deflection. 
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For some relation to give q (pitching velocity), con- 
sider the motion of the airplane as it executes a turn. 
Using vector notation, 


&= pita t+rk 
where @ = angular velocity of coordinate. 


If ¥ = change of direction in horizontal plane from 
original flight path (degree of turn executed), then, 


a= w+ & sin 0} + * cos a 
Therefore, r = (dV/dt) cos ®, g = (dW/dt) sin , and 
dv/dt = rsec ® (9) 
“Eliminating dv /dt, 


g=rtan® (10) 


Using assumptions 1, 2, and 3, and g = + tan ®, the 
longitudinal equations may be written: 


uX, + wX,, — g(w/Uo) = 0 
uZ, + wZ, + 6Zo, = — Ug 
uM, + wM,, + 6,.Mo, = (D — M,)q 
X,, = 0*; therefore w = 0, and so 
uZ, + Zo,6, = —Uoqg 
uM, + Mo,6, = (D — M,)q 
Solving for 6,: 
6, = — (OM,q+ Z,(D — M,)q)/(Zo.M, — MoZ,) (11) 
Integrating Eq. (9) gives the turn as a function of 


time 
‘ 
v= J r sec ®(t)dt 
0 


Knowing W, the acceleration produced during the 
turn can easily be calculated. 
Normal acceleration = U,@ ft./sec.?, where, Q = 


(12) 


adv /dt. 
Gravitational acceleration = 32.2 ft./sec.2 There- 
fore, 
a = V(a,)? + (a,)? (13) 


CALCULATIONS FOR TIME TO EXECUTE A STEADY TURN 
AND THE CONTROL DEFLECTIONS NECESSARY FOR A 
TYPICAL TWIN-ENGINED ATTACK BOMBER 


Dimensions and Characteristics of Airplane Considered 


465 sq.ft. 
61.33 ft. 


Wing area, S 
Wing span, 0 


I 


* X,, = 0 and the assumption X) = 0 are not consistent ex- 
cept n the case for horizontal flight in a straight line at the speed 
U,. The error produced by this inconsistency is small. 
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Aspect ratio of wing, AR = 8.09 
Weight, W . = 19,300 Ibs. 
Taper ratio, o =: 2.66:1 
M.A.C., ¢ = 100.2 in. 
Vertical surface area, S, = 63.28 sq.ft. 
Aspect ratio of vertical 

surface, AR, = 1.28 
Distance from c.g. to 

rudder hinge line, /, =x 307.6 in. 
Distance from c.g. to 

elevator hinge line,/, = 316.5 in. 


Horizontal surface effi- { = 0.60 (power on) 


. ciency, ny | = 0.72 (power off) 
= —0.170 (power on) 
dC/4Cr {= —0.271 (power off) 
Damping factor, K = 5/, 
Horizontal surface area, 
Sy = 101 sqft. 
Aspect ratio of horizon- 
tal surface, AR, = 4,50 
x; = 49.2 ft.? 
} = 59.8 ft.? 
K2 = 104.1 ft.? 
Cow = 0.035 
U, J = 280 m.p.h. (sea level) 
: \ = 300 m.p.h. (15,000 ft.) 
rT.’ J = 0.0304 (sea level) 
” 1 = 0.0395 (15,000 ft.) 
J = 0:87 (sea level) 
| 


Propeller efficiency, 7 0.85 (15,000 ft.) 


Resistance and Rotary Derivatives for Longitudinal Mo- 
tions 





X, = —(2/mU))q’SC, 
Z, = —(2/mU))q'SC, 
1, \? a 
M, = (q'S/mUs)(K)(S,/S (=) : 
@ q~/ 0) (K)(S,/S) K, 1 + (@/"AR, 


The calculation of M, presents some difficulty. If 
moment curves are available for various values of 
T;'(Fig. 1), it is possible to get an approximate value 
of M,,. 

It is known that T,’V* = T¢,'Vo? = const. or, re 
writing, V = Vo(To,’/Tc’)'", where ( )o indicates 
equilibrium trim. 

For power on, Tg,’ = Pn/q’SVo, where Pn = power 
available, g’ = dynamic pressure, S = wing area, and 
Vo = velocity for trim. 

By interpolating a Cy vs. C, curve for T¢,’ and shift- 
ing the C, axis to intercept the curve at C;,, values of 
Cy’ for various values of 7,’ at trim C,; may be ob- 
tained. Since V is a function of T¢’, it is therefore pos 
sible to plot various values of V for corresponding 
values of C,,’, 


uM, = Cy’q’St(g/w) 
M,! = q'St(g/w)(Cy’/u) 
M,,! = (1/C,,)gt(dCy'/dV) 
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TABLE 1 
SrasBitity DERIVATIVES AND ACCELERATIONS FOR UNIT CONTROL MOVEMENT FOR TWIN-MoToRED ATTACK BoMBER 





Vor Zu Mr Loe Nor Noo 


Altitude Ne & Mg My Ly Np Lr Ny 
Sea level* 0 -—0.151 —2.83 0.00201 —7.00 -—0.0694 0.326 -—0.771 —0.579 —41.8 0.467 0.771 0.1775 0.00443 
15,000 ft. t 0 0.109 —2.41 —0.00157 —4.71 —0.0468 0.221 -—0.521 —0.420 -—30.3 —0.339 0.559 0.1289 0.00325 





oc- = O00 } My, based on change in Cy’ at these constant values of Cz, 
+t C- = 0.287 f “* . st ; 


TABLE 2 
TWIN-ENGINED AIRPLANE; SEA LEvEL CONDITIONS 





r sec. ® WV, deg. 6,, deg. 5,, deg. a,, ft./sec.* a,, ft./sec.* 
0 0 11.48 0 32.2 
0.01 0.32 11.99 1.60 
.07 0.64 11. 5.08 
.20 0.77 10. 9.84 
42 0.77 : 15.01 
.14 0.68 5. 26.32 
.27 0.56 ; 38.34 
81 0.43 .16 62.98 
.00 : . 08 85.57 
AT 3 : 99.59 
.81 110.66 
119.35 
121.36 


9 o9 
od 


bt bo bh bo be 


mw De 
Wo 0 ww Ww WO Oo 
NYO KWKHHKHD 


eo 
4 


IWNNNNYNYNHNYNHNNNWNHDWD 


we 
bo 


| 
| 





TABLE 3 
TWIN-ENGINED AIRPLANE; 15,000-Fr. CoNDITIONS 





«t—9? 


@, rad. r r sec, ® W, deg. 6,, deg. 54, deg. a, ft./sec.? a,, ft./sec.* g 

0 0 0 0 0 15.83 32.2 1.00 
0.056 0.00376 .0038 .O1 .85 13.40 : 32.2 1.00 
0.168 0.01152 .0117 .07 yf} 11.28 5. 32. 1.01 
0.310 .02147 .0225 .19 .89 .42 32. 1.05 
0.449 .03093 .0343 .40 . 87 32. 1.10 
0.694 0.04613 .0601 .07 .74 32. 1.29 
0.878 .05572 .0874 13 .58 32 1.56 
1.102 0.06488 . 1484 5.43 41 32 2.20 
1.212 0.06818 .1943 RY £33 32. 2.84 
1.260 0.06938 . 2269 .30 .30 32. 3.26 
1 3 
1 3 
1 


© 


. 289 .07005 .2521 3.16 29 32. 5 
. 308 0.07042 ).2719 38.18 .28 32 
.313 .07049 . 2765 .29 32. 


.85 
3.91 


DdydyONNHNHNHNNHNWN 


eeooorwusn 





where dCy’/d V may be obtained by taking the slope of | where 

the Cy,’ vs. V curve at Vo (Fig. 2). da, _ _ Cow _o. [ly \*Sy dry , f 4 
Dividing M,’ by K,? gives the desired quantity, d(rb/2V) i. Wee | eae (reference 4) 

M,. 


Accelerations for Unit Control Movement 
Resistance Derivatives for Lateral Motion 
Yo, = Cy,.q'(S/m) 


1\ /1 6 \? dC j 4 
L, = (?) (7) ( ) FE oi | (reference 3) " Caf pete q 
T 4 K, d(pb/2V) Mo, = Cuo,9 S(t/m)K,? 
Lo C9’ S(b/m)K,? 


(r =m/p SV) . 
No = Cw9,9'S(b/m)K,? 


b \? d I 'S(b/ 2 

) | so | (reference 3) Nog Ca, S(6/m)K, 
K, d(pb/2V) 

Control Defleciions and Transition Time Necessary for 

Airplane to Maneuver from Straight and Level Flight 


2) Fe | (reference 3) into a Steady Turn 
Re d(pb/2V) ae ; 
An examination of the exponential curve 
: y | 
K,/ |d(pb/2V) 





n n-+m 


a ou — yp~(nt+m)t 
of | (reference 3) =k ( =< i ) 



















Po: 0 
6-084 
Po-.178 
P¢-285 
Po-396 















































Fic. 1. Power on pitching moment coefficient curves. 
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JO \ -lFOR $EA LEVEL |(U.=280 MPH 
Vac av: a8 x|10°5 
08 For 15,000 FT [= 300 MPH) 
dCu/AV= 1.00 X |104 
06 a oes: 
Cu \ 
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Fic. 2. Effect of power on pitching moment coefficient at 
trim lift coefficient. Cz 
15,000 ft. 


trim 
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ANGLE OF BANK & VELOCITY 
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Fic.3. Assumed velocity of roll and angle of bank. 


will reveal that the parameter, k, determines the maxi- 
mum angle of bank, whereas the other two parameters, 
n and m, determine the rolling velocity and, therefore, 
the helix angle pb/2V. 

For a maximum bank of 75° and a helix angle of ap- 
proximately 0.090, the parameters should be k = 2.95, 
m = 3.0, and ” = 1.5. ; 

Using the values of the stability derivatives and ac- 
celerations for unit control deflection for sea level (see 
Table 1), A = 126.5, C, = —0.1102, and C, = 0.1618. 


= ().209 at sea level and 0.287 at 
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Fic. 4. Time to turn and necessary control deflections for twin- 
engined airplane. 
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Fic. 5. Acceleration during turn for twin-engined airplane and 
single-engined airplane. 
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Fic. 6. Time to turn and necessary control deflections for single- 
engined airplane. 


By assuming various values of ¢ and graphically inte- 
grating Eq. (12), the time to enter into steady turning 
flight can be obtained. Substitution of the proper 
values into Eqs. (4) and (5) gives control deflections, 
6, and 6,, necessary to execute this maneuver. Since 
d¥/dt is known from Eq. (9), the acceleration produced 
can be readily calculated (Table 2). 

Fig. 3 shows ® and p plotted against time ¢, while 
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TABLE 4 


STaBILIty DERIVATIVES AND ACCELERATIONS FOR UNIT CONTROL MOVEMENT FOR M 9pIFIeD SINGLS-ENGINED AIRPLANE 





Altitude Xe Zu p My 


—0.00201 —12.37 


Sea level* 0 
15,000 ft.t 0 


—0.109 





* = 
CL 0.588 } M,, based on change in Cy’ at these constant values of Cr. 


tCL = 0.287 


Fig. 4 shows the time to turn along with the necessary 
aileron and rudder deflections. The acceleration pro- 
duced during the maneuver is shown in Fig. 5. 

Eq. (11) gives 6, necessary to hold steady turning 
flight. Since g = r tan ®, 6,,,,,, can be obtained by 
using a large value of time, ¢. For this particular case, 
5, = —8.4°. 

Considering the same values of k, m, and m and using 
the proper stability derivatives and accelerations for 
unit control deflection for 15,000 ft., the time to enter 
into a steady turn and the control deflections necessary 
can be calculated (see Table 3). 

Fig. 4 shows (dashed lines) the time to turn along 
with the necessary aileron and rudder deflections, and 
Fig. 5 shows the acceleration produced during the turn. 

6, = —10.5° in the steady state turn. 

From an examination of Fig. 4 it is apparent that 
there is little difference in the time required to make a 
specified turn assuming ample aileron control. As 
would be expected, however, the control deflections 
must be considerably increased for the high-altitude 
case, which might make it impossible for the pilot to 
obtain the same rolling velocity at the two altitudes, 
meaning a loss in maneuverability as the airplane as- 
cends. 

The elevator deflection necessary to hold steady turn- 
ing flight is slightly increased due to altitude, but the 
acceleration produced remains unchanged. 

The maneuvering characteristics investigated have 
been for a twin-motored attack bomber of late design. 
By assuming a modified single-engined airplane with 
one motor mounted in the fuselage and of the same 
weight and total horsepower as the two motors on the 
twin-engined airplane, the aerodynamic qualities neces- 
sary for the calculations remain unchanged. Only the 
radii of gyration come in as influencing factors. There- 
fore it should be of interest to investigate the maneuver- 
ing qualities of an aerodynamically similar airplane, 
but one with different mass distribution. 

Assume that the top speed remains unchanged and 
the single motor is mounted 6 ft. forward of the two 
motors on the previously discussed airplane. 

The known quantities are: 

(1) Radii of gyration of twin-motored 
K,? = 49.2 ft.4, K,? = 59.8 ft.4, and K,? = 

(2) Weight of engine on twin-motored 
W = 2690 lbs. 

(3) Weight of engine on single-motored 
2W = 5380 lbs. 


airplane: 
104.1 ft.? 
airplane: 


airplane: 


ly Np 
—0.0779 0.576 
—0.00157 — 8.33 —0.0675 0.391 —0.585 


L, N Yor Zo Mee Loa Nor Noe 


1.360 0.1994 0.00669 
0.986 0.1447 0.00490 


—0.467 
—0.339 


—41.8 
—30.3 


—0.579 


—0.420 


—0.865 


(4) Distance of motor from fuselage on twin-motored 
airplane: y = 8.75 ft. 

(5) Distance motor moved forward on single-mo- 
tored airplane: x = 6.00 ft. 

By denoting J; and J; as moment of inertia of single- 
and twin-motored airplane, respectively: J; = mass of 
airplane X Kx? — inertia due to motors. 


Therefore, 
Kx? = I;/mass of airplane 


In a similar manner K,? and K,* can be approximated 
for the single-engined airplane. 

The new values for radii of gyration are: 

27.9 ft.2, K,? = 69.8 ft.2, and K,? = 92.9 ft.? 

Using these new radii of gyration, the stability deriva- 
tives and accelerations for unit control deflection 
can be calculated (see Table 4). 

Again considering the same values of k, m, and m and 
using the proper stability derivatives and accelerations 
from Table 4, the time to maneuver into steady turn- 
ing flight and the control deflections necessary can be 
calculated. Fig. 6 shows the time to turn and the neces- 
sary aileron and rudder deflections for sea level and 
15,000 ft. Fig. 5 shows the acceleration produced. 


5, = —8.4° at sea level in the steady state turn. 
6, = —10.5° at 15,000 ft. in the steady state turn. 


It is of interest to note that for this ‘‘new’’ airplane 
consisting of one engine the time to execute a given 
turn remains substantially the same as that required 
for the original twin-engined one. The initial aileron 
deflection at sea level for this single-engined airplane is 
fairly low, but from Fig. 6 it is seen that the aileron 
position must reach in short time almost that required 
for the first airplane considered. However, at 15,000 
ft. the initial aileron deflection is increased only a few 
degrees and reaches a maximum value after a finite 
time period, which differs from the aileron position- 
time curve for the twin-engined airplane, which shows a 
maximum deflection at zero time. For this modified 
airplane the maximum aileron position necessary at 
15,000 ft. is only a fraction of one degree greater than 
that necessary at sea level, whereas, for the twin-en- 
gined airplane, the altitude case required over 4° addi- 
tional deflection. 

The elevator positions are the same for both airplanes 
at a given altitude and accelerations unchanged with 
regard to altitude and configuration. 

The rudder deflections required are small for all cases. 
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TABLE 5 
MobpIFIED SINGLE-ENGINED AIRPLANE; SEA LEVEL CONDITIONS 





t, sec. ®, rad. r r sec ® W, deg. 
0 0 0 0 0 
0.125 0.056 0.00402 0.0040 0.01 
0.250 0.168 0.01237 0.0125 0.07 
0.375 0.310 0.02306 0.0242 0.21 
0.500 0.449 0.03319 0.0368 0.42 
0.750 0.694 0.04945 0.0644 1.15 
1.000 0.878 0.05970 0.0936 2.29 
1.500 1.102 0.06950 0.1537 5.83 
2.000 1.212 0.07301 0.2081 11.01 
2.500 1.260 0.07429 0.2430 17.48 
3.000 1.289 0.07499 0.2699 24.82 
4.000 1.308 0.07541 0.2912 40.90 

© 1.313 0.07549 0.2961 














6,, deg. é,, deg. 4,, ft./sec.* ay, ft./sec.* ““g 
0 6.50 0 32.2 1.00 
0.26 9.81 1.64 32.2 1.00 
0.52 10.66 5.13 32.2 1.01 
0.63 10.25 9.92 32.2 1.05 
0.63 9.34 15.09 32.2 1.10 
0.57 7.06 26.40 32.2 1.29 
0.49 4.83 38.38 32.2 Rs 
0.39 2.42 63.02 32.2 2. 
0.35 1.15 85.32 32.2 2. 
0.33 0, 52 99.63 32.2 3. 
0.34 0.22 110.66 32.2 3. 
0.32 0.03 119.39 32.2 3. 
0.32 —0.03 121.40 32.2 3. 








TABLE 6 






MobpIFIED SINGLE-ENGINED AIRPLANE; 15,000-Fr. CoNDITIONS 











t, sec. @ rad. r r sec ® W, deg. 
0 0 0 0 0 
0.125 0.056 0.00380 0.00380 0.01 
0.250 0.168 0.01163 0.01179 0.07 
0.375 0.310 0.02162 0.02270 0.19 
0.500 0.449 0.03107 0.03449 0.40 
0.750 0.694 0.04624 0.06020 1.08 
1.000 0.878 0.05580 0.08749 2.14 
1.500 1.102 0.06493 0.14356 5.45 
2.000 1.212 0.06819 0.19434 10.29 
2.500 1.260 0.06939 0.22697 16.32 
3.000 1.289 0.07005 0.25211 23.19 
4.000 1.308 0.07042 0.27189 38.20 

rc) 1.313 0.07049 0.27646 





CONCLUSION 


The effect of altitude (15,000 ft.) on the maneuvering 
characteristics of a twin-motored airplane is apparent. 
There appears to be little difference in the time to make 
a specified turn, but the aileron deflection must be 
considerably increased for the high-altitude case. 

At sea level, the investigation of the maneuvering 
characteristics of the single-engined airplane—i.e., the 
twin-engined airplane considered to be powered by a 
single engine but otherwise unchanged—showed that 
for this modified condition the time to make a given 
turn and the rudder and aileron deflections necessary 
to execute this maneuver differed little from that re- 
quired for the original twin-engined airplane. The 
initial aileron deflection required for the modified air- 
plane is less than that necessary for the twin-engined 
one but quickly builds up to nearly the same maxi- 
mum. Therefore, at sea level it appears that there is 
little difference in the maneuverability of the two air- 
planes. 








é,, deg. 5,, deg. a, ft./sec.* a,, ft./sec.* _ 
0 8.95 0 32.2 1.00 
0.30 11.38 1.67 32.2 1.00 
0.61 10.30 5.19 32.2 1.01 
0.73 9.48 9.99 32.2 1.05 
0.72 8.43 15.16 32.2 2.21 
0.62 6.22 26.49 32.2 1.30 
0.50 4.10 38.50 32.2 1.87 
0.36 2.09 63.17 32.2 2.20 
0.32 0.99 85.51 32.2 2.84 
0.29 0.45 99.87 32.2 3.26 
0.2 0.19 110.93 . 32.2 3.59 
0.29 0.02 119.63 32.2 3.85 
0. 121.64 32.2 3.91 


28 —0.02 





The effect of altitude on the single-engined airplane 
is slight. The initial aileron deflection increases, but 
the maximum deflection required remains the same. 
The rudder deflections and time to make a specified 
turn are only slightly changed. Since altitude has 
such a small effect on this modified single-engined air- 
plane, it certainly should be superior in maneuvering 
qualities to the twin-engined one. 

Accelerations remain the same for all cases consid- 
ered. 
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Fatigue Characteristics of Natural and Resin- 
Impregnated, Compressed, Laminated Woods 


FORREST B. FULLER* AND TURE T. OBERGT 
Wright Field 


SUMMARY 


Because of the increased demand for the production of air- 
craft, natural woods are being used for primary structural mem- 
bers and propeller blades are being fabricated from both natural 
and resin-impregnated, compressed, laminated hard woods. 
Although there are considerable data available on the static 
properties of wood, information about action under repeated 
stresses is quite limited.'~® Therefore, fatigue properties have 
been determined for natural and compressed woods on specimens 
with and without notches. The rotating beam fatigue machine, 
which employs a constant load and has been satisfactorily used 
for testing metals, was also found suitable for testing wood. 
Tests were run to 100,000,000 reversals or cycles, at relatively 
low and high frequencies, a few being continued to 500,000,000 
cycles. Relatively low notch sensitivity was found for the 
natural and compressed woods. The ratio of the fatigue limit in 
bending to the modulus of rupture in static bending was about 0.2 
for circular or round sections. 


MATERIALS 


yo HARD MAPLE and yellow birch woods of 
aircraft grade with straight grain, as well as non- 
impregnated and resin-impregnated, compressed, lami- 
nated hard maple and yellow birch of commercial proc- 
essing, were tested. A description of the materials, 
including specific gravity and moisture content, is 
given in Table 1. The values for moisture content for 
the compressed, resin-impregnated woods are in reality 
percentage loss of weight rather than moisture content. 
The impregnating and bonding materials were thermo- 
setting phenol-formaldehyde resins of the paper film 
and liquid types. Laminated, as used herein, refers 
to that type of construction in which the laminations 
or veneers (thin pieces of wood) are bonded together 
with the grain in all laminations parallel as contrasted 
with plywood in which the grain in alternate lamina- 
tions is parallel but perpendicular to that in adjacent 
laminations. 


TEST PROCEDURE 
Static Tests 


Tensile, compressive and bend specimens were pre- 
pared and tested in accordance with procedures that 
are now incorporated in Federal Specification L-P- 
406.° The bend specimens were rectangular in cross 


Received February 3, 1943. 

* Major, Air Corps, Army Air Forces Matériel Center. 

} Aeronautical Materials Testing Engineer, Army Air Forces 
Matériel Center. 


81 


TABLE 1 
Description of Material 





Specific 
Gravity 


Moisture 
Content, 
Per Cent 


Species Description 





Natural wood—solid (nonlami- 
nated) 

Natural wood, '/g in. lamina- 
tions, noncompressed, nonim- 
pregnated, phenolic resin 
paper bonding 

Compressed wood, '/g in. lami 
nations, 50 per cent, impreg- 
nated, phenolic resin content, 
12 per cent 

Compressed wood, '/,_ in. 
laminations, 50 per cent, 
impregnated, phenolic resin 
content, 23 per cent 

Compressed wood, !/g in. lami- 
nations, 50 per cent, phenolic 
resin paper bonding, non- 
impregnated 

Natural wood—solid (nonlami 
nated) 

Compressed 
laminations, 
35 per cent, 
phenolic resin content, 
per cent 

Compressed wood, !/, in. lami 
nations, 50 per cent, impreg- 
nated, phenolic resin con- 
tent, 23 per cent 

Compressed wood, !/g in. lami- 
nations, nonimpregnated, 
phenolic resin bonding 


8.0 0.70 Hard maple 


3.3 0.68 Hard maple 


Hard maple 


Hard maple 


Hard maple 


Yellow birch 
wood, in 
approximately 

impregnated, 
35 


Yellow birch /is 


Hard maple 


Yellow birch 





* This specimen was the large type having 1 in. diameter in the gauge. 


section with a span approximately 16 times the depth. 
All specimens were tested under comparable laboratory 
conditions of temperature and humidity and with the 
stresses occurring parallel to the grain of the wood. 
Moisture content and specific gravity were determined 
in accordance with standard practice for wood as 
specified in Army-Navy aeronautical specifications on 
the oven dry weight and volume basis." 


Endurance or Fatigue Tests 

Tests were made on ball-bearing rotating-beam 
fatigue machines, similar to the R. R. Moore type, 
Fig. 1. In this machine the specimen is a rotating 
beam of circular cross section supported at the ends and 
loaded with two symmetric loads producing a constant 
bending moment between the loads. The specimen is 
rotated under load, the stresses alternating equally 
between tension and compression. The unit stress in 
the extreme fiber of the specimen is computed by the 
standard bending formula S = MC/I, where 
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Rotating beam fatigue testing machine—wood specimen 
in place. 


Fie. 1. 


S = modulus of rupture or the unit stress in the 
extreme fiber in compression or tension, 
Ibs. per sq.in. 

M = bending moment of the section, in.-Ibs. 

I =moment of inertia of the cross section, in.* 

C = distance from the neutral axis to the extreme 
fiber or the radius at mid-span, in. 


The type of specimens, the type of surface finishing 
and the testing were also in accordance with pro- 
cedures that are now incorporated in Federal Specifi- 
cation L-P-406. Tests were made at frequencies of 
3,450 and 10,600 r.p.m. on notched and unnotched speci- 
mens of the short type, Fig. 2. Stresses for notched 
specimens were based on the root section of the notch. 
The length of the specimens was parallel to the grain 
of the wood. The reduced section of each specimen 
after turning in a lathe was finished successively with 
sandpaper and Nos. 1, 0 and 00 emery paper. 


ROTATING BEAM FATIGUE SPECIMENS 
FOR PLASTIC MATERIALS 





















Fic. 2. 


The tests were conducted in a room in which the 
temperature and humidity conditions did not fluctuate 
appreciably. The weight of the specimens after test- 
ing was not determined, but it is believed that no seri- 
ous change in moisture content occurred during the 
tests. 7 

Impact data were obtained on '/2-in. square speci- 
mens of the standard notched Izod type. with the 
notch transverse to the grain and parallel to the face 
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of the laminations. The specimens were tested in an 
Olsen plastics impact tester of 100-in.lb. capacity. 


RESULTS 


The mechanical properties of the various woods are 
given in Table 2. The data are incomplete for some 
of the woods because of the lack of sufficient material. 

The fatigue or endurance limit is that stress to which 
a material may be subjected for an indefinite number of 
reversals without failure. While the fatigue limit for 
steels in general is obtained at 10,000,000 stress cycles, 
it has been found that a larger number is necessary for 
aluminum alloy. In an investigation on the fatigue 
strength of natural walnut, ash, pine and spruce in the 
rotating-beam type fatigue machine, the knee on the 
S-N (stress-number of cycles) curves was reported be- 
tween 20,000 and 2,000,000 stress reversals—that is, 
an endurance or fatigue limit was determined between 
20,000 and 2,000,000 stress reversals.+* 


TABLE 2 


Average Static Mechanical Properties Parallel to the Grain 











Ultimate Ultimate Bending 
Tensile Compressive Modulus of 

Izod Impact, Strength, Strength, Rupture, 

Test Ft.Lbs. per Lbs. per Lbs. per Lbs. per 
No. In. of Notch Sq.In. Sq.In. Sq.In. 
1 9,700 19,600 
2 14,200 19,000 
3 re i 30,000 24,000 41,000 
4 6.0 33,000 27,000 41,000 
5 41,000 
6 11,200 21,700 
7 3.11 31,300 31,000 
8 8.6 25,000 25,600 32,300 
9 24,400 48,200 










Fatigue data for the natural and compressed woods 
are given in Table 3 and shown in Fig. 3. The data on 
the S-N curve for resin impregnated, compressed, 
laminated yellow birch indicate that the actual fatigue 
limit was not obtained at 500,000,000 stress cycles, 
one specimen failing in the fatigue machine at 435,000,- 
000 cycles at a stress of 10,000 Ibs. per sq.in. The tests 
on the specimens for this particular material were 
stopped at 500,000,000 stress cycles, although this was 
not sufficient to reach the knee on the S-N curve. Re- 
tests were made on these specimens by increasing the 
stress without stopping the machine and the specimens 
run to failure. The number of cycles for each retest 
was compared with that for the specimen tested initi- 
ally at the same stress. If this value for the retested 
specimens was lower than that for the original or initial 
test, the material was probably injured during the initial 
stressing. Only the original data are plotted on the 
S-N curve, Fig. 3, the values for the retest being given 
in Table 4. 

For unnotched material, the fatigue limit for natural 
hard maple wood was determined at 4,000 Ibs. per sq.in. 


















FATIGUE CHARACTERISTICS OF WOODS 








'S-N DIAGRAMS FOR REVERSED BENDING TESTS OF 
IMPREGNATED LOW RESIN CONTENT SO% COMPRESSED 
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STRESS LBS. PER S$Q.IN. x 1000 


NUMBER OF CYCLES OF STRESS FOR RUPTURE - LOG SCALE 


S-N_ DIAGRAMS FOR REVERSED BENDING 
TESTS OF IMPREGNATED HIGH RESIN 
CONTENT, 50% COMPRESSED 
HARD MAPLE, 
NG SYMBOL “TYPE OF SPECIMEN “OPERATING SPEED RPA 
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TABLE 3 


Comparison of Endurance Limit with Static Properties 





Endurance Limit, Rotating Beam, Lbs. per 
Sq.In. “ ven 

Notched specimen 
operating speed, 
——F 3. 
3,450 10,600 


Unnotched specimen 
operating speed, 
r.p.m. 

3,450 10,600 
4,000* 
4,000* 
9,500T 
10,000t 
9,500T 
6,000* 


Z 
9 


3,000* 
8,0007 
9,000¢ 
8,000¢ 


9,500T 
9,500 Tt 


7,500T 
8,000T 


NOP WD 


9,000f 9,000f 8,000f 


oo 


7,500** 


9 10,000* 9,000* 


Ratio of Endurance Limit of Unnotched Material to Average 
Ultimate Static Properties 
Bending 








Circular 
cross section 
0.173 
0.178 
0.197 
0.197 
0.197 
0.234 
0.246 
0.197 
0.176 


Rectangular 
cross section 
204 


Compression 
.412 0. 

. 282 0.210 

.396 0.232 

.302 0.232 

. 232 

.276 
. 290 
. 232 


. 207 


Tension 


. 536 


. 293 
.410 





* Endurance limit based on 100,000,000 cycles. 
+ Endurance limit based on 200,000,000 cycles. 
t Endurance limit based on 500,000,000 cycles, 


** Endurance limit based on-100,000,000 cycles and 1200 r.p.m. 


and for natural yellow birch at 6,000 Ibs. per sq.in., 
although the moisture content and specific gravity were 


about the same. The values for the compressed woods 
were higher with correspondingly higher specific grav- 
ity. No appreciable variation in fatigue strength was 
found between compressed maple containing 23 and 12 
per cent phenolic resin. _In the high-resin material, 


longitudinal cracks in the lamination were noted in 
both the unnotched and notched specimens after 200,- 
000,000: stress cycles, Fig. 4. The characteristics or 
appearances of the fracture in the natural and com- 
pressed woods are different from those of metals. 
In general, they are irregular and jagged throughout 
the cross section, Figs. 6 and 8, whereas in metals 
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TABLE 4 


Fatigue Values of Impregnated, Compressed Yellow Birch Wood 
—Unnotched Specimens 





——Original Tests—-—~ (ae inner Eg 

Stress, Stress, 

Ibs. per Number of Ibs. per Number of 
sq.in. cycles sq.in cycles 
12,000 45,000,000 F* 

12,000 91,100,000 F* 

12,000 32,000,000 F* 

12,000 15,200,000 F* 

10,000 500,100,000 12,000 8,889,000 F * 
9,500 515,200,000 12,000 10,600,000 F* 


9,000 500,350,000 12,000 15,600,100 F * 








* F denotes failure of the specimen in the fatigue machine. 





Fic. 4. Note cracks in lamination in compressed, impreg- 
nated, phenol-formaldehyde, high-resin content, laminated, hard 
maple, notched fatigue specimen. 


there is a smooth surface due to the closing and open- 
ing of the crack during the cyclic stressing, followed 
by a rather small irregular area occurring at the time 
of complete failure. 

For notched material, the fatigue limits were some- 
what lower. However, the impregnated, compressed, 
laminated maple and birch are low in notch sensitivity, 
being relatively insensitive to the effect of notches 
under repeated bending stresses. No data of this type 
were obtained on natural wood. Notch sensitivity is 
the ratio of the fatigue limit of a notched specimen to 
that of an unnotched specimen. The variation in 
fatigue strength with increase in the resin content 
from 12 to 23 per cent for the compressed maple for 
notched specimens was negligible as noted in Table 3. 
The notched impregnated, compressed, laminated yel- 
low birch failed in shear under repeated bending stres- 
ses. This type of failure is shown in natural size in 
Fig. 5 and in an enlarged view in Fig. 6. The failures 
of the other materials are in general normal. 

The results given in Table 3 and Fig. 3 indicate no 
appreciable difference in the fatigue limits at a fre- 
quency of 3,450 and 10,600 r.p.m. for impregnated, 
compressed, laminated wood. No data were obtained 
on natural wood. However, the sloping part of the 
curves in the S-N diagrams, Fig. 3, which represents 
actual failure of specimens in the fatigue machine, 
indicates that the materials have lower fatigue re- 
sistance at stresses above the fatigue limit at the higher 
frequency. 
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Hotched Notched 


lmpregnated, High Resin Content, SOf Compressed 
Yellow Birch 


Dudcal Fatieue sractuce in Laminated Bood Snack 


Bote: Longitudine] shear fatigue feilure in notched specinens 
in yellow birch wod. 


Fic. 5. 


In metallic materials a definite relation has been 
found between the fatigue limit in bending and the 
static ultimate tensile strength. This ratio for ferrous 
metals varies from 0.44 to 0.55 and for nonferrous 
metals from 0.22 to 0.45. For woods, the ratios of the 
bending endurance limit on round specimens to static 
tensile, compressive and bend data are given in Table 3. 
The static bend data that were obtained on rectangular 
specimens have been corrected for circular sections by a 
form factor. The ratios, based on tensile and com- 
pressive data, are quite erratic, but those for the bend- 
ing are more uniform. The natural and impregnated, 
compressed birch had higher ratios than maple. 

A series of fatigue tests at 1,200 r.p.m. was also made 
on large specimens, identified as Test No. 8, approxi- 
mately 1 in. in diameter in the gauge section compared 
to the small specimens from compressed, impregnated 
hard maple with 23 per cent resin content. The mate- 
rial was similar in construction and fabrication to that 
used for the small specimens. The machine, compar- 
able in principle to the small R. R. Moore machine, 





Fic. 6. View of fatigue failures in compressed, impregnated, 
phenol-formaldehyde high-resin content, laminated yellow birch 
jagged in the unnotched and shear in the notched specimen. 



































FATIGUE CHARACTERISTICS OF WOODS 


Rotating beam fatigue testing machine. Capacity— 


approximately 2,000 Ibs. 


Fic. 7. 


Fic. 8. Rotating beam fatigue failure. Large specimen. Com- 
pressed, impregnated hard maple with resin content of 23 per cent. 


Fig. 7, was of larger capacity throughout, could be 
loaded either direct or by a lever system, and was driven 





TABLE 5 


Fatigue Values of Impregnated, Compressed Hard Maple Wood— 
Large Unnotched Specimens (See Fig. 8) 








—Original Tests —— — 
Stress, 

Ibs. per sq.in. 
12,000 
10,000 

8,000 
7,500 
7,500 


Number of cycles 
3,201,000 F* 
4,450,000 F* 

82,500,000 F* 
101,000,000 
105,000,000 





* F denotes failure of the specimen in the fatigue machine. 


though the fatigue limit was lower than that for small 
specimens, the mechanical properties were also lower. 
However, it is of interest to note that the ratios of 
fatigue limit to static bending were comparable, Test 
No. 4 and No. 8, Table 3. A typical failure is shown 
in Fig. 8. The fatigue limit using five specimens was 
determined at 100,000,000 cycles, two specimens being 
stressed at 7,500 Ibs. per sq.in. without failure and one 
at 8,000, one at 10,000, and one at 12,000 Ibs. per sq. in. 
with failures, Table 5. Sufficient tests have not as yet 
been made to obtain a complete fatigue curve; how- 


by a variable speed motor of 600 to 2,100 r.p.m. Al- ever. the data for the five tests are indicative 
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Effect of Cold-Work and Heat-Treatment 
on 24ST Aluminum Alloy 


BRUCE MITCHELL* 


INTRODUCTION 


Sec COMPRESSIVE YIELD STRENGTH of 24ST alumi- 
num alloy extrusions is listed in the ANC-5 as 
38,000 Ibs. per sq.in. This value, generally used for 
design, is the guaranteed minimum for this material, 
but the actual yield strength of some specimens may 
run as high as 56,000 Ibs. per sq.in. Thus 18,000 Ibs. 
per sq.in. additional strength is lost to the designer be- 
cause his ‘‘weakest link”’ has a yield strength of 38,000 
Ibs. per sq.in. 

Aircraft engineers have been testing and studying 
this material for some time in an endeavor to utilize 
this possible increase in strength. The growing list of 
papers published on this subject indicates some gain in 
cold-working for additional strength, an increased 
knowledge of this alloy in the plastic region, and some 
attempts at predicting the compressive yield from the 
tensile strength of test specimens. 

All of the data used in the preparation of this paper 
were obtained from tests made in the laboratory of Con- 
solidated Aircraft Corporation. From these data, a 
more accurate method of calculating the compressive 
yield from the tensile properties will be developed and, 
it is hoped, some additional light will be thrown on the 
general characteristics of 24ST aluminum alloy. 


DESCRIPTION OF TESTS 


During the course of this investigation at Consoli- 
dated more than a hundred tests were made on 54 dif- 
ferent hat-type stringers under a number of different 
conditions oi heat-treatment and cold-work. Forty of 
these were extruded hats 11/2 in. deep with a nominal 
gauge of 0.078 in. The remaining 14 were rolled 
hats 11/2 in. deep made from 24ST Alclad sheet. The 
Alclad hats were made from standard sheets in gauges 
ranging from 0.072 to 0.128 in. Sections of the hats 
5'/, in. long were riveted together, as seen in Fig. 1, 
for the compression tests. 

All compression tests were made in a 200,000-lb. 
Baldwin-Southwark hydraulic testing machine with a 
1,000:1 Huggenberger extensometer mounted on each 
crown of the test specimen. The elongation was taken 
as the average reading of these two extensometers. The 
compression specimens, with their ends ground flat 
and parallel, were mounted on a self-aligning base in 
order to avoid eccentric loading. é 
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Tests of prestretched extruded hat stiffeners compression 
setup. 


Fie. 1. 


All tension tests were made in a 5,000-Ib. Tinius- 
Olsen hydraulic testing machine. Elongation was 
measured with a Riehle extensometer, which can be 
read to '/1o90 in. The tension’ specimens were of 
standard shape with the minimum section 5/s in. 
wide. 

All heat-treating was done in a standard heat-treat- 
ment bath, with a difference in soaking time as the only 
variation of the heat-treatment. This bath is a solu- 
tion of sodium nitrate and potassium nitrate held to a 
constant temperature of 930°F. Quenching was done in 
water at room temperature. After heat-treatment, in 
every case, the specimens were aged 24 hours or more 
before being stretched and three days or more before 
being static tested. 

Stretching was the only method used for applying 
cold-work in this particular group of tests. 

The conventional 0.2 per cent of permanent set was 
used as the criterion in determining the yield strength. 


DISCUSSION 


The mechanical properties of extruded hats differ 
widely, both in the ‘‘as-received” condition and after a 
normal 45-min. heat-treatment. This difference was 
thought to be due to a variation of chemical properties 
until a check of two of the hats revealed almost iden- 
tical quantities of the various constituents. Further 
heat-treatment of these two hats finally equalized 
them after nine hours of soaking time. The weaker 
of the two hats had an ultimate tensile strength of 
65,000 Ibs. per sq.in. after the first 45-min. heat-treat- 
ment. It still tested the same after nine hours. The 
other hat had an ultimate tensile strength of 77,000 









24ST ALUMINUM 


TABLE 1 
Mechanical Properties of 24ST Extruded Hat Stretched Succes- 
sive Amounts 


Stretch, 


Ur Yr Ye 
71,000 46,700 45,700 
69,000 51,000 45,300 
71,500 53,500 45,700 
72,500 55,000 46,600 
73,500 56,000 47,000 
72,800 57,700 48,300 
75,250 61,500 49,300 
77,500 68,000 54,500 
80,600 73,200 57,800 





Ibs. per sq.in. after the first 45 min. and 65,000 Ibs. per 
sq.in. at the end of nine hours. Since the strength of 
the hats can be equalized by a prolonged heat-treat- 
ment, a uniform base line can be established for deter- 
mining the amount of cold-work in any specimen. The 
establishment of a base line is a start, but the type of 
cold-work must also be known before the mechanical 
properties can be properly analyzed. 

Cold-work increases the tensile yield strength of the 
material quite rapidly, but the ultimate unit strength 
is increased only in direct proportion to the change in 
cross-sectional area of the section. The tensile yield 
strength will approach the ultimate as the amount of 
work is increased until they are almost equal, then cold- 
working is stopped just short of the breaking point. 
Thus, the difference between the ultimate and the yield 
is an indication of the amount of work in the material. 
This difference averaged about 23,000 in all test speci- 
mens after heat-treatment. However, some of the 
cold-work may remain after heat-treatment, as evi- 
denced by the nine-hour soaking time required in the 
case of the one specimen previously mentioned. The 
cold-work that remains in the material after a normal 
heat-treatment will be referred to here as “residual” 
cold-work in order to distinguish it from the cold-work 
that reduces the difference between the ultimate and 
the yield. In other words, residual cold-work is the 
amount of strength an extrusion has after a normal 
heat-treatment above what it would have had if given 
a prolonged heat-treatment. 


SYMBOLS 


= tensile ultimate lbs. per sq.in. 
tensile yield Ibs. per sq.in. 
compressive yield lbs. per sq.in. 
23,000 — (Up — Yr) 
Ye when X = 0 
Y, when X = 0 


DEVELOPMENT OF FORMULA 


The compressive yield strength of these extrusions 
after a prolonged heat-treatment was found to be 


ALLOY 


/2 


Compressive yield curve of 24S-T extruded hat stretched 
successive amounts. 


Fic. 2. 


O 4 8 12 1€& 
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Fic. 3. Tensile yield curve of 24S-T extruded hat stretched 
successive amounts. 
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Fic. 4. Residual cold-work curve of 24S-T extrusions. 


42,000 Ibs. per sq.in. The total compressive yield 
strength at any time, regardless of its previous treat- 
ment, will therefore be 42,000 plus residual cold-work 
plus compressive gain due to other types of cold- 
work. 

Since stretching was the only method used to cold- 
work these test specimens, it is the only type of cold- 
work combined with residual cold-work in these speci- 
mens. 

Table 1 is a list of the mechanical properties of a 
single extruded hat that was first heat-treated and then 
stretched successive amounts. The compressive yield 
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TABLE 2 
Mechanical Properties of 24ST Extruded Hats 
” Error 
: bes Ye Ye Per Xo oa7 * <o17 
Ur Yr Ua > Ve xX 19,000 (Calc.) (Act.) Cent Up Ur 

65,400 41,800 23,600 —600 19 43,500 43,700 + 0.5 —0.01 
63,600 42,500 21,100 1,900 190 42,750 44,000 + 2.8 0.03 
65,400 42,700 22,700 300 5 44,300 43,200 — 2.6 0 
63,100 42,800 20,300 2,700 384 42,300 43,800 + 3.4 0.04 
61,500 44,300 17,200 5,800 1,770 41,450 44,000 + 5.8 0.09 
62,100 44,800 17,300 5,700 1,710 41,800 43,300 + 3.5 0.09 
61,400 44,400 17,000 6,000 1,895 41,400 42,100 + 1.7 0:10 
68,000 47,300 20,700 2,300 279 45,300 46,500 + 2.6 0.03 
63,500 48,700 14,800 8,200 3,540 43,400 43,200 — 0.5 0.13 
64,600 49,200 15,400 7,600 3,040 43,800 42,900 — 2.1 0.12 
63,300 49,000 14,300 8,700 3,980 43,500 42,500 — 2.4 0.14 
77,000 50,000 27,000 —4,000 842 45,750 46,700 + 2.0 —0.05 
66,500 49,500 17,000 6,000 1,895 44,500 46,500 + 4.3 0.09 
77,100 49,600 27,500 —4,500 1,065 45,250 47,500 + 4.7 —0.06 
70,200 50,000 20,200 2,800 412 46,550 48,200 + 3.4 0.04 
73,300 51,000 22,300 700 26 48,900 50,800 + 3.9 0.01 
71,600 51,000 20,600 2,400 303 47,450 47,200 — 0.5 0.03 
72,150 52,800 19,350 3,650 700 47,650 47,200 — 1.0 0.05 
69,600 52,900 16,700 6,300 2,090 46,400 45,300 — 2.4 0.09 
71,400 52,800 18,600 4,460 1,020 47,200 47,200 0 0.06 
74,900 52,500 22,400 600 19 49,900 50,000 + 0.2 0.01 
74,700 52,800 21,900 1,100 64 49,600 51,200 + 3.1 0.01 
74,400 53,300 21,100 1,900 190 49,250 49,000 — 0.5 0.03 
75,600 52,800 22,800 200 2 50,500 50,300 — 0.4 0 
74,750 54,000 20,750 2,250 266 49,350 49,000 — 0.7 0.03 
74,000 54,000 20,000 3,000 474 48,800 48,000 im Bod 0.04 
75,100 54,800 20,300 2,700 384 49,500 50,300 + 1.6 0.04 
76,100 54,200 21,900 1,100 64 50,450 50,500 + 0.1 ere 0.01 
65,400 54,600 10,800 12,200 7,830 47,150 43,500 — 8.4 0.19 
62,750 55,500 7,250 15,750 13,050 49,300 43,600 —13.1 0.25 
62,500 56,400 6,100 16,900 15,070 50,650 44,800 —13.1 0.27 mt 
68,000 56,300 11,700 11,300 6,725 48,000 47,400 —- 1.3 hate 0.17 
66,500 56,800 9,700 13,300 9,320 48,850 47,400 — 3.1 0.20 
67,200 57,700 9,500 13,500 9,600 49,450 48,300 — 2.4 0.20 
67,200 57,900 9,300 13,700 9,880 49,650 49,000 —- 1.3 0.20 
70,000 59,200 10,800 12,200 7,840 49,900 52,300 + 4.6 0.18 
66,500 59,700 6,800 16,200 13,800 52,100 55,500 + 6.1 0.24 
70,100 60,500 9,600 13,400 9,450 51,100 52,000 + 1.7 0.19 
74,900 60,700 14,200 8,800 4,075 50,500 50,300 — 0.4 0. 
74,100 61,500 12,600 10,400 5,690 51,000 49,600 — 2.8 <7: 0.1 
70,500 61,900 8,600 14,400 10,900 52,300 50,000 — 4.6 0.20 
69,000 62,300 6,700 16,300 14,000 56,050 49,000 —14.4 0.24 Hee 
76,900 62,700 14,200 8,800 4,075 51,700 54,100 + 4.4 pes 0.11 
71,700 65,300 6,400 16,600 14,500 55,750 52,300 — 6.6 0.23 
71,600 65,600 6,000 17,000 15,200 56,200 50,400 —10.4 0.24 esa 
76,500 65,500 11,000 12,000 7,580 53,600 53,000 — 1.1 0.11 
77,000 66,000 11,000 12,000 7,580 53,900 56,000 + 3.7 ace 0.11 
74,000 66,700 7,300 15,700 12,980 56,000 50,000 —12.0 0.21 

strength is plotted against ““X’’ in Fig. 2. The resultant Y, = Yr, + 1.7X (2) 

curve is a parabola that may be written as: and for negative values of “X” 

Ye = (X*/19,000) + Yo (1) Vr = Yz, — 1.7X (3) 


Yo is the residual cold-work plus the basic strength, 
and X?/19,000 is the stretch type of cold-work. 

The Y; values of Table 1 plotted against ‘“X”’ pro- 
duce the curve in Fig 3 where for positive values of 
“x 





It can be seen from Fig. 2 that the difference be- 
tween Y, and Y¢ is very smali for values of “X”’ up 
to +1,000. Since VY» is the compressive yield when 
“X” is zero and since Yz, is the tensile yield at the 
same point, the connecting link between Y) and Y, 
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TABLE 3 
Mechanical Properties of 24ST Alclad Hats 








: Error 
r ‘ See Yo Ye Per * onz * <ouz 
Ur Y, Us — Ye x 19,000 (Calc.) (Act.) Cent Up Ur 

65,600 47,300 18,300 4,700 1,160 43,750 44,200 +1.0 0.07 
67,750 51,000 16,750 6,250 2,060 45,300 44,000 —3.0 0.09 
68,100 52,000 16,100 6,900 2,510 45,700 45,000 —1.6 0.10 
66,600 52,500 14,100 8,900 4,070 43,500 44.600 2.5 in 0.13 
66,900 56,200 10,700 12,300 7,960 48,150 46,500 —3.5 0.18 oon 

67,900 53,300 14,600 8,4C0 3,710 46,150 46,300 +0.3 0.12 
66,400 53,300 13,100 9,900 5,160 46,050 45,300 —1.7 0.15 
68,000 53,800 14,200 8,800 4,075 46,350 46,200 —0.3 0.13 
68,100 54,400 13,700 9,300 4,550 46,700 45,500 —2.6 0.14 
67,800 54,400 13,400 9,600 4,850 46,700 47,400 +1.5 0.14 
68,300 56,000 12,300 10,700 6,025 47,700 47,800 +0.2 0.16 
68,750 56,000 12,750 10,250 5,530 47,650 49,600 +3.9 0.15 
68,600 56,300 12,300 10,700 6,025 47,900 48,500 +1.2 0.16 
69.000 56,400 12,600 10,400 5,700 47,950 46,700 —2.7 0.15 

can be obtained by plotting Y) against Y7,. Choosing LIMITATIONS OF FORMULA 


values of “‘X’’ less than 1,100 from Table 2 and plotting 
the corresponding Y7, (calculated from Eqs. (2) and (3)) 
values against Yo, the curve in Fig. 4is obtained. The 
equation for this curve is: 


Y = 0.6Y7, + 19,000 (4) 


Combining Eqs. (1), (2), and (4), and Eqs. (1), (3), and 
(4) 
Y- = (X2/19,000) + 1.02X + 0.6¥, + 19,000 (5) 


The sign of the term “1.02X” in the above equation is 
negative when ‘‘X”’ is positive and positive when ‘““X” 
is negative. The product of this term is, therefore, 
always negative and Eq. (5) can be written as: 


Yo = (X*/19,000) — 1.02X + 0.6Y7 + 19,000 (6) 


if the numerical value only of ‘““X”’ is used. 


PROOF OF FORMULA 


_ Tables 2 and 3 list the actual tensile and compres- 
sive strengths, together with the calculated compressive 
strength, of all of those hats on which both tension and 
compression tests were run. Table 2 has a few cases 
where the error in calculating the compressive yield 
is quite high. However, only those with both a low 
ultimate and high value for X have this large error. 
On this basis a limiting factor can be set up by using the 
ratio of X/Ur. This ratio was set at X/Ur = 0.17 asa 
safe limit. This allows a maximum error of 5.8 per 
cent in this series of tests. As a further check on the 
probable cause of error, a 20-ft. length of extruded hat 
was tested for variation in tensile strength. It was 
first heat-treated and stretched 1'/2 per cent. Then 
ten tensile samples were cut from sections throughout 
its length. The tensile ultimate of these ten speci- 
mens varied +4 per cent and the tensile yield +3 per 
cent from the average. 


Because it is an expression of residual cold-work 
and the stretch type of cold-work, this formula will not 
apply to precompressed specimens. However, it can 
be used for determining the original compressive 
strength of a specimen that has failed in compression. 
The high compressive stress, applied just before failure, 
on such a specimen will increase the compressive yield 
strength but will not appreciably affect the tensile prop- 
erties. 

Since the thickest hat section tested in the develop- 
ment of this formula was 0.128 in., and the thinnest 
0.072 in., it should not be used outside of these limits, 
and because of the variation of tensile properties in a 
single stringer, it cannot be expected to give results 
closer than +4 percent. It should not be relied upon 
if Xy, is greater than 0.17. 


APPLICATION OF FORMULA 


In order to take advantage of the full strength of ex- 
truded hats in aircraft construction, the following pro- 
cedure is suggested: All extrusions intended for use in 
stringer combinations would be numbered and weighed 
separately and tensile samples cut from each. The 
tensile samples would be tested and the compressive 
strength of each extrusion calculated by means of Eq. 
6. (The calculated value should be reduced 4 per cent 
to allow for material variations in the stringer.) The 
average area of each extrusion can be calculated from 
its weight and length. Then the actual compressive 
strength in pounds would be figured and the extrusions 
sorted according to their respective strengths. 

The designer, knowing the actual strength of the 
stringers, would place the strongest stringers in the 
region of highest stress. For example, the shear lag 
across a wing puts a higher load on the stringers nearest 
the spars. These stringers would therefore be taken 
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from the high-strength group and the weaker ones 
would be placed farther from the spars. 


CONCLUSIONS 


24ST aluminum alloy extrusions have a potential 
compressive yield strength of 56,000 lbs. per sq.in. but, 
because of variation in mechanical properties, have a 
design allowable strength of only 38,000 Ibs. per sq.in. 
This design allowable can be increased up to 30 per 
cent if a tensile test is made of each extrusion and the 
true compressive strength of the extrusion determined 
by Eq. 6. Then, by weeding out the weakest sections 
or placing them at points of lowest stress, the maxi- 
mum compressive yield strength of the plate stringer 
combination can be realized. 


Letter to 


Dear Sir: 

An inspection of Eqs. (1) to (4) of my paper entitled ‘‘Pro- 
portioning a Canard for Longitudinal Stability and Safety 
Against Stall,’ published in the December, 1942, issue of the 
JOURNAL OF THE AERONAUTICAL SCIENCES, leads to a further 
conclusion of particular interest which was not listed at the end 
of that paper. 

In order to do away with the limitation of tail length, and 
therefore, of range of allowable c.g. locations, which is one of the 
major disadvantages of the canard, the tail surface should be 
designed in such a way that its lift coefficient would be in- 
dependent of the angle of attack of the airplane. In this case 
the derivative of the term containing the tail length and tail 
surface characteristics would be zero, and these items would be 
eliminated from Eqs. (3) and (4). Then the static longitudinal 
stability of the airplane would be independent of the tail length 
and the most rearward allowable position of the c.g. for stability 
would be as far back as it would be in the tailless airplane. 

This could be attained by attaching the horizontal tail sur- 
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the Editor 


faces to the body of the airplane by means of an articulation that 
would leave the tail surface free to rotate about a transversal 
axis. Then the horizontal tail surface, if so designed as to be 
longitudinally autostable, would always remain at its angle of 
attack of equilibrium, regardless of the angle of attack of the rest 
of the airplane. The pilot, through the control stick and tabs or 
other auxiliary surfaces attached to the tail, would control the angle 
of attack of equilibrium of the latter, not its incidence with re- 
spect to the rest of the airplane. 

This type of canard would still be automatically safe against 
stall if h, (Eq. (5)) was assumed as the most rearward c.g. loca- 
tion. However, since such a horizontal tail surface would be 
fully controllable at all angles of attack of the airplane, ip could 
be assumed as the most rearward allowable c.g. location; this 
would provide a further increase of the range of allowable c.g. 
displacements. 

Josepa V. Foa 
Chief Engineer 
American Aero-Marine Industries, Inc. 














Graphical Analysis of Delay of Response in 
Air-Speed Indicators 


KALMAN J. DeJUHASZ*+ 
The Pennsylvania State College 


SUMMARY 


The delay of response of the pneumatic system of the air-speed 
indicator is analyzed by a graphic method and illustrated by an 
example of an actual installation, with various simplifying as- 
sumptions and for various modifications of the system. The 
results are compared with those previously obtained by mathe- 
matical analysis. The results of mathematical analysis agree 
well with those obtained by graphic analysis for those pneu- 
matic systems in which the inertia of the moving air column is 
small in comparison with the resistance to flow; in other cases 
important deviations have been found. The design data for 
minimum delay have been determined. Applicability of the 
method for any given law of impressed pressure change is illus- 
trated by examples. 


—_ LIFT OF AN AIRPLANE is dependent on the air 
speed; therefore it is of prime importance that the 
pilot should be informed at all times of the accurate 
value of the speed. For this purpose there is used an 
air-speed indicator,! which is based on the principle 
that the difference of pressure Ap between the dynamic 
pressure at a pitot head and the static pressure of the 
surrounding atmosphere is proportional to the square 
of the air speed v, according to the equation 


Ap = pv,”/2 


where p is the density of the air. 

The pressure difference is measured by a sensitive 
pressure gauge, usually of the aneroid capsule type, 
which is mounted on the instrument board in the fuse- 
lage. The pitot element is placed in the undisturbed 
air stream at a distance of several feet from the pressure 
gauge. The two are connected with a tubing which, in 
actual installations, is about 30 to 50 ft. long. A damp- 
ing orifice is sometimes incorporated in the system in 
order to damp out vibrations. These elements con- 
stitute a pneumatic system, illustrated in Fig. 1, com- 
prising volumes and resistances to flow. Because of 
the presence of these elements, any change of pressure 
occurring at the pitot head will not be communicated 
to the pressure gauge instantaneously but with a certain 
delay. Another source of delay is the inertia of the 
moving parts of the pressure gauge itself, but this is 


Presented at the Radio and Instruments Session, Eleventh 
Annual Meeting, I.Ae.S., New York, January 25-29, 1943. 

* Professor of Engineering Research, School of Engineering. 

{ The author expresses his appreciation to Drs. W. G. Brom- 
bacher and C. S. Draper for their comments and suggestions and 
to Mr. A. Yorgiadis for the careful drawings and for checking 
the calculations. 
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small in comparison with the delay caused by the pneu- 
matic system and will not be considered in the present 
treatment. 

The purpose of this paper is to (1) determine the 
delay in terms of the design data of the pneumatic 
system and (2) find the dimensions of the pneumatic 
system whereby the minimum delay is attained. 


STATEMENT OF THE PROBLEM 


In the present treatment the altitude and tempera- 
ture, and hence the static pressure and density, are 
assumed to be constant. In order to simplify the prob- 
lem it is assumed that the air speed changes suddenly, 
as in a gust, involving an instantaneous change of dy- 
namic pressure at the pitot head from an initial value 
p: to a new value p2 (Fig. 2). Then the pressure p in 
the pressure gauge will not assume the new dynamic 
pressure instantaneously but will approach it gradually 
as some function of time p f(t); hence during a 
certain time it will not indicate the air speed cor- 











rectly. Itis our purpose to determine the function p = 
f(t). 
P 
m, 
D B f — ——= 
Fae 
Pood 
L ;* > 
Pa 
/ 
/ 
A | t 
on =~ 
Fic. 1 Fic. 2 
Fic. 1. Pneumatic system of air-speed indicator. P, pitot 
head; D, damping orifice; LZ, tubing; A, aneroid capsule. 
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“. 


Response of air-speed indicator to sudden pressure 
change. Dynamic pressure increases suddenly from value ~; 
to value ~,. The air-speed indicator indicates a pressure p = 
f(t), which may be a periodic or aperiodic function of time. 


Fic. 


This problem has been treated in the engineering 
literature experimentally and by mathematical analysis, 
both for related aircraft instruments**~® and specifi- 
cally for air-speed meters.**® In the usual mathemati- 
cal analysis, however, rather drastic simplifying as- 
sumptions are involved. In the present paper a 
graphic method will generally be used which permits a 
closer approximation to the actual conditions and pos- 
sesses the advantage of simplicity and clarity of pres- 
entation of the phenomena involved. Similar graphic 
methods have been applied to other problems of tran- 
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sient phenomena in elastic systems, but these have not 
been applied to problems of aircraft instruments so far 
as it could be ascertained. 


Basic RELATIONSHIPS AND ASSUMPTIONS 


(a) The change of state of air is assumed to be 
adiabatic, for which (p/p) = p/po; dp/po = dp/kpo; 

= 1.4; p) = 10,000 kg./m.? In the present case the 
pressure change is assumed to be small in comparison 
with the absolute value of pressure; therefore the refer- 
ence pressure fp is assumed to be constant. The value 
kPyp = 14,000 kg./m.? is the “modulus of volume elas- 
ticity”’ of the air. 

(b) A container is characterized by its “‘capacity”’ 
C, which is the volume rate at which the medium enters 
into it or leaves it per unit of pressure change in the con- 
tainer and which can be expressed by the equation: 


._ 4Vp) _ pdV | Vp 
pdp pdp  pdp 


For a rigid container dV/dp = 0; therefore the capacity 
is: 





Crigia = V/xpo 
For an elastic container, the volume of which is a linear 
function of the pressure: 


Vesastic ~ Vo + (dV/dp)p mas Vo + Ap 


the capacity can be expressed as 


(Vo/xpo) + A = 


Caastic fas 


a (Vi ot+ Akpo) - 
The term Veguivatent = Vo + Axpois the volume of a rigid 
container having the same capacity as the elastic con- 
tainer and Will be termed “equivalent rigid volume.”’ 
(c) The pressure drop due to an orifice or pipe is a 
function of the flow velocity Ap = f(v), which is a simple 
proportionality in the case of laminar flow. In the 
present case the flow velocity is low and the diameters 
are small; in other words, the Reynolds number is low; 
therefore laminar flow is assumed to exist for which the 
pressure drop is expressed by the Hagen-Poiseuille 
formula: = (32 L/d*)v; L =.length; d = diame- 
ter; v = flow velocity; » = absolute viscosity of air. 
The volume flow per second, g, can be expressed: 


q = vF = (d*F/32 nL)) Ap = Ap/R; 
R = 32 nL/d?F = 
Symbols with index Z will refer to the tubing; with 


index a, to the aneroid capsule; and with index D, to the 
damping capillary orifice. 


KPo ( Veqetvetent ) 


resistance of pipe 


MATHEMATICAL ANALYSIS 


It is assumed that all the volumes, capacities, and 
resistances in the system are combined into one total 
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volume, total capacity, and total resistance, for which- 
V=V,+ V,; C = V/«eh; R= R, + Rp. Then 


for a q volume flow per second it can be written: 
Ap = qR; qdt = Cdp 


Broadly, then, the mathematical analysis treats the 
problem in the simple form: A compact container C 
with a negligible length dimension, Fig. 3a, communi- 
cates with the outside atmosphere through an orifice 
having the resistance R; the pressure of the outside 
atmosphere changes instantaneously from p; to pr»; 
the pressure p in the container is to be determined as a 
function of time p = f(t). Thus, the inertia of the 
flowing air column is neglected, and it is assumed that 
at any given instant the pressure at all points in the 
pneumatic system is the same. It can be written: 
Cdp = [(pe al p)/R\dt; hence, dp ‘(pe —_ p) == dt/CR. 
After integration the exponential relationship is ob- 
tained: 


bz — p = (pz — pie /*° 


which is depicted in Fig. 3b. 


Rh fe ee Pa 
a ie 
re & 
/ | Qay 
p 


(a) Substitute pneumatic system as treated by 
(b) Change 








Fic. 3. 
mathematical analysis; C, capacity; R, resistance. 
of pressure in substitute pneumatic system. 


The product RC = + has the dimension of time and is 
termed the ‘‘time constant”’ or “characteristic time’’ of 
the system. It characterizes the rapidity of response 
of the system to a change of pressure, and its significance 
is: that after the elapse of time 7 the error, Ap, will be 
1/e = 0.368 times its initial value. (After the elapse of 
time 27, the error is 13 per cent; after time 37, it is 5 
per cent; after time 4r, it is 1.8 per cent; and after 
time 57, it is 0.67 per cent.) 

It will be noted that this formula is analogous to 
that describing the electrical phenomenon of charging 
and discharging a condenser through an ohmic resist- 
ance: For a sudden change of the terminal voltage 
from £, to EF, it can be written: 


k,-E= (Ee a Ey)e~ “8° 


whereby p2 — p is analogous with the voltage differ- 
ence, C with the capitance, and R with the resistance 
in the circuit, Fig. 4. 

It is evident that the mathematical analysis repre- 
sents a rather rough approximation to the actual condi- 
tions, because it assumes in the system an infinitely 
great velocity of propagation of the pressure change 
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Fic. 4. Electrical analogy of pneumatic system; EF, E2, im- 
pressed voltage; R, ohmic resistance; C, condenser; EF, voltage 
across condenser. 


and neglects the effect of the momentum of the air 
column. Mathematical difficulties prevent taking 
these effects into consideration by mathematical 
analysis, but it can be done with comparative ease by 
graphic analysis. 


GRAPHIC ANALYSIS 


The equalization of pressure in air-speed indicators 
belongs in the group of transient phenomena in elastic 
columns. These have been treated by graphic analy- 
sis in previous publications’® where also the full deriva- 
tion is given. For this reason at this place only the 
main conclusions, as applied to the present problem, 
will be summarized. 

The graphic analysis deals with ‘“‘disturbances,” 
which term denotes a change of velocity or of pres- 
sure, impressed on a fluid column at some point in the 
system as a function of time. In the present case the 
instantaneous change of pressure at the pitot head, 
from ~; to pz is a disturbance that is representable 
in a t-p diagram as a step function. 

In graphic analysis the elements of the system and 
the phenomena are represented as functions of two of 
the following variables: pressure p, velocity v, time /, 
and location x. The following relations hold good: 

(1) In a system a disturbance, originated at one 
point x» at a certain instant of time f, is propagated 
along the system with a constant velocity a termed the 
“velocity of propagation’’ or ‘‘acoustic velocity.” 
Therefore for the location x of the disturbance at an 
instant ¢ of time it can be written: 


xX — X Ax 


=— = +q = tan (+9) 
i-~ (*9 





This relationship can be represented in a ¢-x chart by 
a straight line forming an angle g with the /-axis. The 
sign is positive or negative depending on whether the 
disturbance travels forward (increasing x with increas- 
ing time) or backward (decreasing x with increasing 
time). For air the mean value of the velocity of 
propagation is: a = 338 m./sec. 

(2) For a finite length of pipe the value L/a = T 
is an important property denoting the time interval 
necessary for a disturbance to travel from one end of the 
pipe to the other. 

(3) In a frictionless pipe a change of pressure is 
always accompanied by a proportionate change of 
velocity and vice versa, i.e., 
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Ap/ Av = const. = +=xpo/a 


It is convenient to use, instead of the linear velocity 

v, the volume velocity of flow g = Fv, with which it 
can be written: 

Ap “) A KPo 


p 
= — = const. = + — 
Aq F Av Fa 


= tan (+a) 


This relationship can be represented in a g-p chart by 
straight lines forming an angle a with the g-axis. The 
expression xpo/ Fa may be termed the “‘wave constant,” 
and its value for air is: 


kpo/Fa = 14,000/338F = (41.5/F) kg. sec./m.5 


(4) The resistance of an orifice or a pipe can be 
expressed by the pressure drop p. — p necessary to 
produce a certain rate of flow gq, i.e., 


(pb: — p)/q = R = tanp 


which can be represented in the g-p chart by a straight 
line forming an angle p with the g-axis and passing 
through the g = 0, p = p2 point. For turbulent flow 
the function is quadratic and can be represented by a 
parabola. The function needs not be stated in mathe- 
matical form but can be represented by a curve ex- 
pressing the flow-velocity, pressure relationship, which 
may be an experimentally determined one. 


APPLICATION OF GRAPHIC ANALYSIS TO THE AIR-SPEED 
INDICATOR 


In possession of the above relationships an air- 
speed meter will be analyzed under various assump- 
tions and modifications. The data of the air and of the 
air-speed indicator® (installation in the German Ju 52 
airplane) are: 

Air 

Modulus of elasticity = «cpp = (1.4 X 104/m.?) kg. 

Acoustic velocity = a = 338 m./sec. 

Wave factor = xpo/a = (41.5/m.*) kg. sec. 

Kinematic viscosity = » = (1.83 XK 10-*/m.?) kg. 

sec. 

Density = p = (0.123/m.*) kg. sec.? 

Aneroid capsule 

Volume [at po = 

m.? 
Volume increase [at Ap = (200/m.*) kg.] = AV = 
0.7 X 10-* m.? 
AV/ Ap = A = 35 X 10-" m.5/kg. 


(104/m.? kg.] = Vo = 11 X 10-* 


Equivalent rigid volume = Veguivatent = (11 X 10-* + 
1.4 X 104 X A) = 60 X 10-* m.? 

Equivalent capacity > C, = V equivalent /kpo = 43 X 
10- m.°/mg. 

Wave factor = «po/aF = (xpo/Vequivatent) (L/a) = 
0.69 X L 


Tube 
Length = ZL, = 10m. 





d 





Diameter = d, = 6 X 107? m. 

Cross-sectional area = F, = d,? 7/4 = 28.3 = 10~° 
m.? 

Volume = V, = L,F, = 283 X 10-*m.? 

Capacity = C, = V,/xkbo = 202 X 10-m.5/kg. 

Resistance = R, = 32nL,/d,? F, = (0.575 X 108/ 
m.5) kg. sec. 

Interval of traverse = L,/a = 0.0296 sec. 


Damping capillary 
Length = L, = 0.2 m. 
Diameter = d, = 0.90 X 107% m. 
Resistance = Rp, = (22.7 X 10°/m.°) kg. sec. 


For this system the calculated time constant is: 


r = (Rp + R,) (C, + C,) = 23.27 X 245 XK 10-4 = 
0.57 sec. 
and the pressure function is: 
pr — p = (bx — pile” 

(a) In the graphic analysis, as a first approximation, 
it will be assumed that the volume of the capsule is 
added to the tube in the form of an additional length, 
so that the total volume of the assumed tube is equal 
to that of the actual tube plus the volume of the cap- 
sule. Furthermore, the resistance of the pipe is added 
to that of the damping orifice. In other words, the 
approximation will be similar to that used in the mathe- 
matical analysis, with the refinement, however, that 
the finite length of the system and the inertia of the 
air column are taken into consideration. 

The length of the substitute tube will be L’ = (283 + 
60)/28.3 = 12.1 m. and its time interval L’/a = 0.0358 
sec. The inertia of the air column is represented by 
the angle a for which 


tan (#a)= Ap/Aq = kpo/Fa = (kpo/V)(L/a) = 
T/C = (1.47 X 10°/m.*) kg. sec. 


The resistance of the system comprises that of both the 
damping orifice and the pipe L and is characterized 
by the angle p for which: 


tan p = (p2 — p)/q = R = (23.27 X 10/m.5) kg. sec. 


These values of a and p are represented in the qg-p 
chart, and the construction is carried out as shown in 
Fig. 5. A pressure difference Ap = 100 kg./m.? is 
assumed corresponding to an air-speed increase from 
160 km. per hour to 216 km. per hour. 

On this Fig. 5 are shown: (a) the ¢-x chart on which 
every point represents a certain location in the system 
as defined by the x-ordinate and an instant of time as 
represented by the abscissa. By lines having the slope 
tan (+g), running from the orifice end to the closed end, 
this coordinate system is divided into triangular zones; 
the condition of the velocity and pressure in each zone 
is defined by a similarly numbered point in the q-p 
chart. (b) The g-p chart contains the R = tan p 
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Fic. 5. Graphic analysis of Ju 52 air-speed indicator installa- 
tion, assuming an added tubing length, so that total tubing 
volume equals the system volume. (a) t-x diagram; (b) q-p 
diagram; (c) history diagram of pressure at pitot head (step 
function suddenly rising from p; to 2}, and at capsule 
(pressure rising in steps of 2L/a duration, with steps of decreasing 
magnitude); calculated exponential pressure curve is also 
drawn in. 


line representing the combined resistances and the 
+tan a lines representing the flow conditions in the 
tube. After the elapse of successive L/a intervals the 
disturbance wave arrives alternately at the orifice end 
(characterized by the R = tan p line) and at the closed 
end (i.e., capsule end, characterized by g = 0). These 
successive conditions are defined by points, 1, 2, 3,.... 
in the g-p chart and are valid for the zones marked 
1, 2, 3,....in the ¢-x chart. (c) The history of pres- 
sure chart for the p’tot end of the tube (characterized 
by the step function rising from p; to f2 instantane- 
ously) and the capsule end of the tube (characterized 
by a stepwise rising pressure line, having steps of de- 
creasing magnitude as the time progresses). 

The calculated exponential pressure curve is also 
drawn in this diagram and shows a close argument with 
the graphically obtained stepped line. 

The sum of the steps up to a certain time can also be 
calculated as follows: After the elapse of 2L/a = 2T 
time the pressure at the closed end of the tube will in- 
crease by 2 Ap, for which 

tana tana = 


Api = [(b2 — pr) — Api] aa a 








therefore the ‘“‘error’’ will be 


b2 — Per = [(1 — m)/(1 + m)] (pb: — prs) 
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Similarly, after the elapse of 47 time the error will be 
bz — Per = [(1 — m)/(1 + m)]? (b2 — pri) 

and after the elapse of ¢ = 2n7 time the error will be 
bz — Por = [(1 — m)/(1 + m)]"(b2 — pri) 


It can be shown that as the length of the container is 
reduced (keeping the volume constant), i.e., the inter- 
val T is reduced, then in the limiting case the stepped 
function will assume the form of an exponential func- 
tion. Namely, substituting T = ¢/(2n) in the above 
equation: 


_, _ f1— @/2 CR) aT re 
a” F + (t/2 CR) (1/n) (bs — Ps) 


is obtained, and in the limiting case, lim L 





= 0; lim 





IT = 0; lim n = o; therefore 
2 — (t/2 Dale) f nite 
1 + (¢/2 CR)(1/n) 
hence 
pr — p, = (be = pie" 


which is the same expression as that obtained by the 
mathematical analysis. 

(b) A closer approximation to the actual conditions 
can be realized by treating the capsule as a distinct 
element (i.e., not to absorb its volume in the tube 
volume, as it was done under (a)). In Fig. 6 it is as- 
sumed that the capsule is represented by a tube having 
the same length as the connecting tubing and having a 
cross-sectional area such that the volume is the same as 
the capsule volume. The wave factor of the capsule 


will be: 
L, = 10m. 
A , ; 
tan B = | hg (kpo/ V) = (6.9 X 10°/m5.) kg. sec. 
q 


The construction is shown for the ¢-x chart, the q-p 
chart, and for the history chart for the pressure, for 
both the pitot end and the capsule end. It is seen 
that the stepwise rising pressure line exhibits a some- 
what different character from that obtained for as- 
sumption (a), but the calculated exponential pressure 
function still closely coincides with the graphic result. 

(c) A still closer approximation to the actual con- 
ditions can be realized by treating the capsule as a tube 
shorter than the connecting tubing. In Fig. 7 it is 
assumed that the capsule length is one-third of the tube 
length, and correspondingly the cross section is greater 
than in the former case. The wave factor will be: 


L. = 3.33 m. 


Ap 
Aq 


a 


tan B = = (2.3 X 10°/m*.) kg. sec. 


The construction is shown in Fig. 7a for the t-x 
chart, Fig. 7b, the g-p chart, and Fig. 7c for the history 
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Fic. 6. Graphic analysis of Ju 52 air-speed indicator installa- 
tion, assuming a capsule length equal to the tube length. (a) 
t-x diagram; (b) g-p diagram; (c) history diagram of pressure 
at pitot head (step function suddenly rising from p; to 2), 
and at capsule (pressure rising and falling in steps of 2L/a 
duration); calculated exponential pressure curve is also drawn 
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Fic. 7. Graphic analysis of Ju 52 air-speed indicator installa- 


tion, assuming a capsule length equal to one-third of the tube 
length. (a) ¢-x diagram; (b) g-p diagram; (c) history diagram 
of pressure at pitot head (step function suddenly rising from ; 
to p2), and at capsule (pressure rising and falling in steps of 2L/3a 
duration); calculated exponential pressure curve is also drawn in. 





; 
{] 
f 
i 





96 JOURNAL OF THE AERONAUTICAL SCIENCES—MARCH, 


chart for the pressure, for both the pitot end and the 
capsule end. It is seen that the stepwise changing 
pressure line exhibits a different character than in the 
formerly treated assumptions (a) and (b), but the 
calculated exponential pressure function still passes 
through the steps of the graphically obtained stepped 
line and can be regarded as the faired, mean value of the 
latter. 

(d) Actually, the length dimension of the capsule 
is a much smaller fraction of the tube length than it has 
been assumed under (c). With the expenditure of 
added labor and attention to the graphic construction 
but without any inherent difficulty, the actual con- 
dition could be realized with any desired approximation 
by the procedure shown in the previous examples. 
This would be, however, only of academic interest be- 
cause even the foregoing examples show that the mathe- 
matical analysis gives a satifactory approximation for a 
system with the design data as given. Nevertheless, 
the graphic construction is of interest, since it gives a 
closer approximation of the phenomenon and, so to 
speak, gives an insight into the microstructure of the 
phenomenon which the mathematical analysis is in- 
capable of providing. 

(e) In the example treated under (a), (b), and (c) 
the relative magnitude of the resistance and inertia 
was such that 


tan a = Ya 
tan p av 


It is of interest to investigate a case of such characteris- 
tics that 


= [(L/a)/(RC)] <1 





tan a/tan p = (T/RC) > 1 


For this purpose the formerly assumed system will be 
treated with the modification that the damping capil- 
lary is removed, so that the only resistance to flow will 
be that caused by the tube itself, i.e., 


tan p = R = (0.575 X 10°/m.°) kg. sec. 


tan a/tan p = 1.47/0.575 = 2.56 
In this case the calculated time constant will be: 
= RC = 0.575 X 245 X 10-* = 0.014 sec. 


The graphic construction is showy in Fig. 8, giving the 
t-x and q-p charts, and also the history chart for the 
pressure, for both the pitot end and the capsule end. 
In the latter the calculated exponential pressure curve 
is also included. It is seen that in this case the actual 
phenomenon will not be an aperiodically damped func- 
tion of time as represented by the exponential pressure 
curve but a periodically damped function as represented 
by the graphically obtained stepped line exhibiting 
steps above and below the true pressure pf, line, the 
steps having decreasing magnitude as the time pro- 
gresses. 

(f) Particular interest attaches to the case in which: 
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tana =tanp; L/a = T= RC 


and which is shown in graphic representation in Fig. 9. 
It is readily seen that in this case complete equalization 
of pressure will occur at the end of the first L/a interval. 
Such a system possesses the least conceivable delay of 
response, which is the time necessary for the ‘“‘new’’ 
air pressure to be propagated to the closed end (i.e., 
capsule end) of the tube. The phenomenon in this 
case is analogous to a critically damped vibration. 
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Fic. 8. Graphic analysis of Ju 52 air-speed indicator installa- 
tion, assuming that the damping orifice is removed and assuming 
an added tube length, so that the total tubing volume equals 
the system volume. (a) t-x diagram; (b) g-p diagram; (c) 
history diagram of pressure at pitot head (step function 
suddenly rising from p; to p2), and at capsule (pressure fluctuating 
above and below the true pressure line 2 with decreasing magni- 
tude as time progresses); calculated exponential pressure curve 
is also drawn in. 
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Fic. 9. Graphic analysis of Ju 52 air-speed indicator in- 
stallation, assuming a wave factor equal to the resistance. (a) 
t-x diagram; (b) g-p diagram; (c) history diagram of pressure 
at pitot head (step function suddenly rising from p,; ,to po), 
and at capsule (step function suddenly rising from p; to p2); 
calculated exponential pressure curve is also drawn in. 


(g) If the volume of the capsule is negligible in 
comparison with the tube volume and, if there is no 
damping resistance, then this optimum condition, 
i.e., least delay, will be obtained when the tube length 
and diameter are in the following relationship: 


L = aRC = a(32 n L/d?F)(FL/kpe) 
d? = 32 n(a/xpo)L 
d=12X10* VL 











AIR-SPEED INDICATORS 97 


For example, for a length of 10 m. the optimum diame- 
ter is 3.8 mm. 

It has to be borne in mind that the above results are 
based on assumptions that may not always be fulfilled 
practically. A streamline flow has been assumed, fol- 
lowing the Hagen-Poiseuille formula, while actually the 
resistance is distributed along the whole length of the 
tube. Therefore, the above results are subject to 
experimental verification and correction. The salient 
conclusion, however, appears to be justified—that for 
minimum delay of response: (1) the volume of the 
pressure gauge should be as small as possible; (2) no 
concentrated resistance, in the form of damping capil- 
lary orifice, should be incorporated in the pneumatic 
system; (3) the tubing diameter should be in the neigh- 
borhood of that stated by the formula d = 1.2 X 
10-* \/LZ (in which the dimensions are stated in meters). 


RESPONSE TO A GRADUAL CHANGE OF AIR SPEED 


In the foregoing analysis a suddenly changing air 
speed has been assumed. With the aid of graphic 
analysis it is possible to determine the response to other 
than instantaneous changes, defined as a function of 
time. As an example, Fig. 10 illustrates the applica- 
tion of the method to a linear change of dynamic pres- 
sure, occurring in 2L/a time, in the pneumatic system 
treated in Fig. 5. It is seen that in this case the steps 
of the response function are not horizontal but slanting, 
forming a continuous broken line with decreasing slope 
toward the time-axis as time progresses. 

With added labor but without any inherent difficulty, 
any other kind of pressure change, such as described by 
a curvilinear history function, can be analyzed with the 
graphic method described. 

The method may be used for the solution of problems 
of other related aircraft instruments, such as altimeters 
and rate-of-climb indicators. Furthermore, it may 
applied to the analysis of engine pressure indicators 
and of the pressure fluctuation phenomena between two 
interconnected spaces, such as the main combustion 
chamber and the air cell of a Diesel engine. In the 
latter cases, however, the pressure difference is great; 
hence, the flow will not be laminar nor can the wave 
factor be regarded as a constant; therefore the pro- 
cedure will be less simple than in the above-treated 
cases. 
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Fic. 10. Graphic analysis of Ju 52 air-speed indicator in- 
stallation under the same assumptions as treated in Fig. 5, 
except that the impressed pressure change is assumed to be 
linear, rising from p; to p2 during 2L/a time. (a) t-x diagram; 
(b) g-p diagram; (c) history diagram of pressure at pitot 
head, and at capsule (pressure rising from ~; in a function of 
continuous slanting broken straight line, with the slope de- 
creasing as time progresses). 
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Wind-Tunnel Cooling 
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SUMMARY 

The problem of cooling wind tunnels has been analyzed. It 
has been found that the turning vanes and the throat of the 
tunnel are the best places for using surface cooling. 

As the operating speeds of tunnels increase, it becomes more 
and more difficult to cool the tunnel surfaces, and it becomes 
necessary to install a radiator in the tunnel. The detrimental 
effect of the use of a cooling radiator is small when the iunnel 
cross section is expanded ahead of the radiator. The percentage 
reduction in the energy ratio of the tunnel is given by the ratio of 
twice the adiabatic temperature rise, when the air in the radiator 
is brought to rest, to the temperature difference available for 
cooling, all multiplied by the ratio of the density of the air in the 
radiator to the density of the air in the throat of the tun- 


nel. 
The length-diameter ratio of the air passages of the cooler is 


determined by the cooling capacity required. The actual size 
of the air passages is limited, on the one hand, by the space avail- 
able and, on the other, by the magnitude of the allowable power 
expenditure pushing the coolant through the radiator. The size 
of the radiator whose coolant power expenditure is equal to the 
power expended in the rest of the coolant circuit will, in general, 
be found to be about the best practical compromise. 


INTRODUCTION 


rr AN EFFCRT to obtain large-scale, high-velocity data, 
the trend in modern wind-tunnel design has been to- 
ward larger tunnels operating at lower than atmospheric 
pressure. But as the size and power requirements of 
wind tunnels increase and as the operating pressures 
decrease, the problem of cooling becomes difficult. An 
examination of the literature reveals a noticeable lack of 
information on how to handle this difficult cooling 
problem. Furthermore, there is no information avail- 
able on the relative merits of various practical solu- 
tions. 

This paper attempts to supply some of this informa- 
tion by considering in detail the two most feasible 
methods of cooling wind tunnels: (1) cooling by heat 
transfer through the tunnel walls and/or the corner 
vanes; and (2) cooling by heat transfer to a radiator 
installed in the tunnel. 


ANALYSIS 


The following is an analysis of the available informa- 
tion applicable to wind-tunnel cooling. 

The equation for heat transfer through a laminar 
boundary layer in a compressible fluid is! 


(d/dt)| JegT + (u?/2)] = (kO*/Oy*)| JepT +( ou?/2)] (1) 
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where 
J = mechanical equivalent of heat 
Cp = specific heat of fluid at constant pressure 
T = temperature 
u = velocity 
k = thermometric conductivity 
y = distance perpendicular to surface 
o = Prandtl number 
t = time 


In the case of heat transfer to a flat plate at a con- 
stant wall temperature, Eckert and Drewitz* found that 
when o = 1, the heat transfer per unit surface per unit 
time is given by 
H/S = ue,(OT/dy)o = ule(Ti — Tw) + (u:?/2J)] X 

(Ou/Oy)o/u1 (2) 


where the new symbols are defined as 


H = heat dissipation 

S = surface 

m = coefficient of absolute viscosity 

(OT /Oy)o = temperature gradient at surface 

T; = temperature outside the boundary layer 
Tw = temperature of surface 

My = velocity outside the boundary layer 


(Ou/Oy)o = velocity gradient at surface 


Eq. (2) should apply very closely for heat transfer 
from air (¢ = 0.73) since the velocity and temperature 
fields in the boundary layer are very nearly the same 
when o = 0.73 and when o = 1. This solution is 
based on the Blasius velocity distribution in the bound- 
ary layer and therefore assumes no pressure drop across 
the plate. It may, however, be regarded as an approxi- 
mate solution to the general heat-transfer equation 
when the pressure gradient along the plate is small. 
This approximate solution will now be put into a very 
convenient form for cooling calculations. 

The pressure gradient at a given section of a wind 
tunnel is given by 


Ap/L = u(0u/dy)o(4/D) 


where D is the diameter of the tunnel at the section 


considered. 
Substituting this relation into Eq. (2), Eq. (3) is 
derived. 


H = {[(T1 + m?)/2I¢] — Tw}(4D?/4) Ap(c,/m) (3) 


Assuming the average velocity at the section is 
approximated by 1, 
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H = {((Ta + tq?)/2Jq|] — Tw}cp>M Ap/pgu,? (4) 


where 
Pq = the average air density at the section 
ua = the average velocity at the section 


T, = the average air temperature at the section 


Eq. (4) is similar in form to the usual heat balance 
equation® encountered in dealing with heat transfer 
from an incompressible fluid to a wall at constant tem- 
perature. 


H _ (T, = Tw)GM(1 _ e~ 4?/pa%a*) 
or at low values of Ap/p,u,” 
H = (Tq — Tw)¢pM Ap/patta’ (5) 


When dealing with heat transfer from a compressible 
fluid to a wall at constant temperature, the usual heat 
balance equation (Eq. (5)) applies, provided the effec- 
tive temperature of the fluid is taken as the actual fluid 
temperature plus the adiabatic temperature rise, u,”/ 
2Jcy, which occurs when the fluid is brought to rest. 
Because Reynolds analogy between the mechanism of 
heat transfer and friction applies to both laminar and 
turbulent air flow, Eqs. (4) and (5) are applicable in 
both cases. 


SURFACE COOLING 


Consider a closed return-circuit wind tunnel with 
small angles of expansion and contraction throughout. 
The flow through such a tunnel will be similar to turbu- 
lent flow in along pipe. The frictional pressure drop is 
therefore approximated by the relation 


ApD/2pqugL = 0.049(pattaD/p)~°-? (6) 


which has been given by McAdams‘ as representing the 
frictional pressure drop in the latter case. The mass 
flow of air through the tunnel being constant, Eq. (6) 
can be put into the more convenient form 


Ap/patta? = (0.098L/D-*)(4M/mp)~%-? (7) 


where the quantity in parenthesis is constant for all 
sections of the tunnel. Eq. (7) shows that the friction, 
Ap/ pata’, per unit length is inversely proportional to 
D°.8, Because the heat dissipation is roughly propor- 
tional to Ap/ pau,” (Eq. 4), the greatest heat dissipation 
per unit length of tunnel will occur at the tunnel throat. 
Because the turning vanes of the tunnel form, in general, 
passages of smaller hydraulic diameter than the throat 
of the tunnel, the vanes offer the most effective surfaces 
for cooling. 

In all well-designed wind tunnels, the working section 
occurs immediately after a rapid contraction in the 
tunnel cross-sectional area. Because of the resultant 
thinning of the boundary layer, the friction is more 
intense in the throat of wind tunnels than in long pipes. 
In other words, the throat of a tunnel is an even more 
efficient place at which to cool than the above analysis 


indicates. Similarly, because the passages between the 
turning vanes of the tunnel are analogous to the en- 
trance length of a pipe,® they will have more friction 
and more cooling capacity than indicated above. 

It has now been established that the working section 
and the turning vanes offer the best surfaces for cooling 
wind tunnels. The next question that arises is how 
much cooling surface will be needed. 

When the cooling at the working section is sufficient 
to carry off the heat generated there, Eq. (8) applies. 


QAp/J = [(Ta + AT,ad) — Ty] (GoM Ap/patta*) (8) 


where Q is the volume flow of air. From Eq. (8) it is 
easily shown that the difference between the air tem- 
perature and the wall temperature must be equal to the 
adiabatic temperature rise of the air when brought to 
rest if sufficient cooling is to be obtained. 


2AT ad = (T, + AT,ad) — Tw (9) 


Eq. (9) has not been further simplified because the 
temperature, 7, + AT,ad, is related to the temperatute 
at which test models will be operating and is therefore 
limited by structural considerations. 

Eq. (9) shows that as the speed at the test section of 
wind tunnels increases it becomes increasingly difficult 
to carry off by surface cooling the heat generated there. 
It is evident that, despite the reduction in the life of the 
turning vanes, they must be used for cooling in all large 
tunnels using surface cooling. The heat dissipation to 
be expected from the vanes can be estimated conserva- 
tively by the use of Eqs. (4) and (6) or, more exactly, by , 
the use of Eq. (4) and the data for the inlet length of 
pipes given by Langhaar.® 

Quite often even these two sets of surfaces—the 
working section and the turning vanes—will not be 
sufficient to handle all the cooling of the tunnel. In 
such cases the next most efficient surfaces for cooling— 
those where the tunnel has small cross-sectional areas— 
will have to be used. 


RADIATOR COOLING 


A wind tunnel can be cooled by means of a radiator in- 
stalled in the tunnel instead of by surface cooling. In 
order to keep the power expended pushing air through 
the radiator and the turbulence of the tunnel at a 
minimum, the radiator should be carefully designed so 
as to have only frictional pressure losses. How much 
bad effect on the performance of a wind tunnel such a 
well-designed radiator would have is easily evalu- 
ated. 

The performance of a wind tunnel is usually described 
by its energy ratio 


E.R. = (Mu;*7/2)/power input (10) 


where the subscript ¢ refers to average air conditions in 
the throat of the tunnel. But all the power put into 
the tunnel eventually appears in the form of heat which 
must be carried off by the coolant system. For ade- 
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quate cooling, therefore, the heat balance of the wind 
tunnel is expressed by 


Mu ?/2J/E.R.= | [T,+(u,?/2Icy)] — Tw} (GM Ap/pru,’) 
or 


((T, + AT ad)/(T, + T,ad)|""—*[2AT,ad + 
(T, + AT,ad)— Ty] = E.R.[Ap/(pm?/2)] (11) 


where the subscript 7 refers to average air conditions in 
the radiator. 

Eq. (11) states that the percentage reduction in the 
energy ratio of a wind tunnel due to the use of a radiator 
is equal to the ratio of twice the adiabatic temperature 
rise of the air in the radiator, when brought to rest, to 
the temperature difference available for cooling, all 
multiplied by the ratio of the density of the air in the 
radiator to the density of the air in the throat of the 
tunnel. Because the temperature, 7, + AT,ad, is 


limited by the allowable temperature of the model and 


the temperature Jy is limited by the coolant system, 
the percentage reduction in the energy ratio of the tun- 
nel depends completely on the air speed in the radiator. 
In other words the reduction in the energy ratio of the 
tunnel can be made as small as desired by expanding 
the cross-sectional area of the tunnel ahead of the 
radiator. 

The adiabatic temperature rise that occurs when the 
air flowing through the radiator is brought to rest is 
easily related to the adiabatic temperature rise equiva- 
lent of the air ahead of the radiator by the geometry of 
the cooler. Thus, in general, the open area of the 
radiator will be about half that of the tunnel cross sec- 
tion in which it is placed. It follows that the velocity 
of the air in the radiator will be twice that ahead of the 
radiator or the adiabatic temperature rise equivalent 
will be four times that ahead of the radiator. Eq. (11) 
can therefore be expressed in the more convenient form, 


E.R.Ap AER. be + re - 
pu?/2 E.R. T, + AT,ad 
SAT ad 
(T, + AT;,ad) — Ty 


where A7,ad equals the adiabatic temperature rise 
equivalent of the air ahead of the radiator. 

The loss in wind-tunnel performance due to the use 
of a radiator for cooling has now been evaluated. It 
has also been shown that in order to have adequate 








(12) 





cooling, a definite value for Ap/p,u,* must be attained 
(Eq. 4). The value of Ap/p,u,? defines the geometry 
of the air passages of the cooler, i.e., the length-diameter 
ratio. The actual radiator depth, however, is a com- 
promise between two factors: the space available and 
the magnitude of the power to be expended pushing the 
coolant through the radiator. The latter decreases in- 
versely as the cube of the depth of the radiator. The 
optimum radiator depth will, in general, be approxi- 
mated by the depth of the radiator whose coolant power 
expenditure is equal to the power expended in the rest 
of the coolant circuit. 


CONCLUSIONS 


The problem of cooling wind tunnels has been ana- 
lyzed. It has been shown that as the operating speeds 
of tunnels increase (1) it becomes more and more 
difficult to cool the walls of the tunnel, and (2) in the 
absence of a radiator, the cooling of the turning vanes 
becomes mandatory. It has also been shown that the 
reduction in the energy ratio of the tunnel due to the 
use of a radiator can be made as small as desired by 
expanding the tunnel cross section ahead of the radia- 
tor. 

In any plan for reducing the effect of the radiator in- 
stallation on the performance of the wind tunnel by the 
use of a combination turning-vane-radiator design, the 
nonturbulent character of the flow will tend to be com- 
promised. Furthermore, at its best the turning-vane- 
radiator is only equivalent, not superior, to an installa- 
tion in which vane cooling plus radiator cooling is used. 
Because of the relative simplicity of construction of the 
latter, the turning-vane-radiator installation should not 
be given serious consideration. 
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The Influence of Sweep 


on the Spanwise 


Lift Distribution of Wings 
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SUMMARY 


A method is developed for determining the twist of a swept- 
back wing with prescribed spanwise lift distribution. Explicit 
formulas are derived for the case that the lift is given at eight 
points of the semispan. 


INTRODUCTION 


oe METHODS for computing the spanwise lift 
distribution of a monoplane wing with arbitrary 
chord distribution have been worked out.! However, 
the lift distribution furnished by these methods agrees 
well with experimental results only if the line containing 
the centers of pressure of the wing sections (approxi- 
mately the quarter chord line) is straight and at a right 
angle to the direction of flight. Sweep back of the wing 
increases the load at the wing tips. This effect has an 
unfavorable influence on the stalling characteristic of 
the wing unless it is compensated for by giving the wing 
a suitable twist. 

The problems connected with the spanwise lift dis- 
tribution of a swept-back wing are of twofold nature: 
Either the shape of a wing is given, and the lift distri- 
bution has to be found, or it may be asked what twist 
has to be given to a swept-back wing in order to get a 
certain lift distribution. The first problem is dealt with 
in papers by Crean? and by Weinig.* The present 
paper is devoted to the second problem, which may be 
formulated more precisely as follows: A straight wing 
with a certain lift distribution may be given sweep 
back by sliding each section back; then find the angle 
of attack at which the swept-back wing sections have to 
be placed in order to preserve the original lift distribu- 
tion. For this new angle of attack a formula is found 
here, which expresses the angle in terms of the given 
lift distribution. In deriving this formula the present 
paper follows closely the procedure of Crean’; but 
whereas Crean’s paper restricts itself to small angles of 
sweep back, no such restriction has been made here. 
Weinig? likewise assumes small sweep back and, fur- 
thermore, considers only wings with rectangular lift 
distribution. 

The coordinate system is taken in the usual way: 
x-axis in the direction of flight, y-axis to starboard, 
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z-axis downward. The wing is assumed to be sym- 
metric with respect to the x-axis. The locus of the cen- 
ters of pressure of the wing sections—approximately at 
one-quarter of the chord for wings without flaps—may 
be called the pressure line. Only wings with pressure 
lines consisting of two straight line segments in the 
x-y plane are considered here. The angle between a 
pressure line segment and the y-axis shall be called the 
angle of sweep A, taken positive for sweep back (see 
Fig. 1). 








Fic. 1. Diagram of a swept-back wing. BCD = pressure line. 


The theory of the lifting line can also be applied 
to the swept-back wing, the lifting line coinciding with 
the pressure line. A lifting line element of length ds 
and vorticity [ along a pressure line swept back through 
an angle A (see Fig. 1) will experience a lift of 


pI'V cos A ds 


V being the velocity in the negative x-direction, and p 
being the density of the air. Since dy = ds cos A, the 
lift of a wing section bounded by two planes perpendic- 
ular to the y-axis is 


dL = pI'V dy 
The vorticity [ is given by the following relation: 
T = agcVao/2 


where dp and ¢ are section characteristics that are un- 
changed when the wing section is slid back; ap is the 
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Fic. 2. A vortex of the swept-back wing. 
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effective angle of attack. The vorticity and likewise the wing. Thus the following result is found: By sweeping 
lift distribution will therefore remain unaltered as long a wing back, the lift distribution will be preserved if the 
as the effective angle ap is not changed. This will be angle of attack of each wing section is changed by an 
the case if angle a*, which is given by 


a — (w/V) = a, — (w,/V) * 


a* = a, — a = (w, — w)/V 
where a denotes the geometric angle of attack and w the 
down-wash velocity for the straight wing, while a, and This additional angle of attack will, of course, vary 


w, are the values of these quantities for the swept-back along the span. 


THE DOWN-WASH VELOCITY 


The vortex system of the swept-back wing is a superposition of vortices of the type shown in Fig. 2. The down- 
wash velocity induced at a point P of the pressure line by the whole vortex system is the sum of the contribu- 
tions that each individual vortex gives. A single vortex ABCDE of vorticity — AT induces at P the velocity 





i 2? + cos ABP , cos PBC + cos BCP , cos PDE + |(- ar) 
An yt y’ PF y= y’ 
It can be seen from Fig. 2 that cos PDE = — sin A; but that is the case only if y’ < y, as in the figure; for 


, the absolute value of y — y’, these two cases can be 





y’ > y, however, cos PDE = + sin A. By using |y — y’ 
combined into one expression: 


cos PDE/(y — y’) = — sin A/|y — y’| 
and Aw, becomes 
1 1 — y’) sin A 1 + y’ cos 2A 
au, = of (1 - o— 7) ) + ee ( Lonnie — cos 24) _ 
4niy + y Vy? + 2yy’ cos 2A + y”? 2y° sin A Vy? + 2yy’ cos 2A + y”? 
sin A 1 1 1 cos 2A sin A 1 
+ + ;|(- ar) - 2 .. aoe - ad . + 
y-yl y-y 4nlLy+y’ 2y’sinA ly—y'] yy 
Vy? + 2yy’ cos 2A + ad (- ar) 
2y’ sin A(y + y’) 





























Since the velocity induced by a straight wing of vorticity — AT is 
Aw = (1/4r){[1/(y + »’)] + [L/(y — »’)]}(— 41) 
the additional down-wash velocity due to sweep back is found to be 


Y in J Vy? + 2yy’ cos 2 2 
al _ cos2A — sind 4 y? + 2yy" cos 2A + y |(-ar)- 




















Aw, — Aw = : | | : : 
4a 2y’ sin A ly — y’| 2y’ sin A(y + y’) 
1 1 I 1 Vy? + 2yy’ + y"? — Ayy’u? 
z[ -so+«(- i) + . + yy FY a ]an 
4m 2uy y psi 2y’u(y + y’) 


where u was written for sin A. 





The term ~/(y? + 2yy’ + y’? — 4yy’u?)/2y’u(y + y’) can be expanded into a series: 


fa el Ses 
2yt (y+y’)? . 29m of “No +7? 


where C, =[1- 1- 3, . . . (2n — 3)]/[2- 4-6, .. . 2m]. In this series each term is less than half the preceding one if 
u? < 1/20rA S 45°. The numerical computations will show that the terms for m = 2 and m = 3 are already 
small; therefore, the terms for n > 3 shall be neglected, so that Aw, — Aw can be written as follows: 


Zaz! Day8ae/2 
aw, — aw =2[ u(L- 1 2) _ yp iy _ yy P| (ar) 
4 yy byl owty)? (y + y’)4 (y + y’)® 


The Additional Angle of Attack 


The vortex system shall be assumed to consist of N vortices. Let yon + 1,m = 0, ..., N — 1, be the span 
coordinates of the points where these vortices branch off the pressure line. Then the down-wash velocity shall be 
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computed at NV — 1 points ya, k = 1,..., N — 1, which lie between the branching-off points of the vortices. If 


the vorticity for the points 2, is given, the vorticity of the vortex branching off at yen, + , is 
— AT = (yom) — T'(yem + 2) 


where it is assumed that I'(0) is also given and that [(b/2) = 0, 6 being the span. Thus the additional down- 
wash velocity due to sweep back induced at y2, by the whole vortex system can be written as a sum: 


1 | ( 1 1 
W—- w= — i + se 
4r Vox lyom oS Yox| 


u3 


Vom + 1 


(Yom +1 + 5) . 





‘ 9 
2Vom + 1° Vox? 


— n° (Vom +1 + a || Pom) =. T'(yom i‘ »| 


Dividing this expression by V, the additional angle of attack a* is obtained. Introducing dimensionless co-, 
ordinates » = y/(b/2) and dimensionless circulation coefficients y = I'/Vb and putting y(nm) = Ym, the addi- 
tional angle of attack a* for the point 2, is found to be 


> 
Yom +1° Vox 


(Yom + 1+ Yar)! 





N-1 
a* (nox) = Dd) (matt — Dmx? — Cet) (Ym — Ym + 1) 


m=0 


where 








1 1 1 N2m + 1 
Omk Ss = — | ra 9 
Qe Linx  |mem+1 — Mexl (mom + 1 + 2x)? 


Dm = (1/27) [nom +2" nox/(N2m t+i1+ nox) *| 
Cm = (1/7) [nom + 1° n2x7/(N2m 41+ nox)® 


This, in turn, can be simplified by putting 


Sox = ox; Fink = Gmzr — Im—-1,% 
Zon = box; Sm = baz — by - lk 
hox = Cox; Im = Cmk —~ Cm— 1,5," = i. sey N-1 


Thus, a*(n,) becomes 
N-1 
a* (nox) - Zz (fF matt - Lmx® — Nnxt®)¥m 
0 


m= 


or 
a*(n,) = Fy sin A — G, sin? A — H, sin® A (1) 
where 
N-1 N-1 N-1 
F, = De SmunYm Gy = i ZmkYm H, = 2d himzYm (2) 
The structure of Eqs. (1) and (2) shows that in order to find a* at the N —' 1 points ny, k = 1,..., N — 1, 


it is necessary to compute only once the coefficients finz, Zmx, and Am,. Then the values F,, G,, and H; for any 
arbitrary lift distribution are obtained by inserting the circulation coefficients y,, of that lift distribution into Eqs. 
(2), and a*—in radians—is given by Eq. (1) in terms of the angle of sweep back A. 


NUMERICAL WORK 





The computation of the coefficients fn, Zmxy and Am, 
was carried out for the case that the vorticity is pre- 
scribed at the eight points of the semispan whose co- 
ordinates are nom = sin (mxr/16),m = 0,...,7. The 
coordinates sin [(2m + 1)2/32] were chosen as 72m + 1. 
The coefficients finz, Zmx, ANd Ay,, m = 0,...,7; k = 1, 
..., 7, are tabulated in Table i, so that a* can readily 













be obtained by means of Eqs. (1) and (2) if the eight 
circulation coefficients yo, . . ., yz are assigned. For a 
set of values 7», which were kindly given by The Glenn 
L. Martin Company, the coefficients F, G, and H are 
tabulated in Table 2, and the angle a*—in degrees— 
for angles of sweep back A = 15° and 30° is plotted in 
Fig. 3. 
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TABLE 1 
k—1 2 3 4 5 6 7 
m Smx 
{ 
0 —1.0054 —0.2107 —0.09788 —0.06028 —0.04364 —0.03538 —0.03142 
1 —0.04691 —1.1978 —0.2823 —0.14290 —0.09575 —0.07768 —0.06539 
2 1.1232 —0.07217 —1.2954 —0.3009 —0.15541 —0.10494 —0.08893 
3 0.2358 1.1669 —0.12825 —1.5158 —0.3684 —0.2011 —0.15015 
4 0.10247 0.2299 1.3045 —0.2263 —1.9156 —0.5063 —0.3093 
5 0.05382 0.09330 0.2301 1.5052 —0.4323 —2.7384 —0.8427 
6 0.02907 0.04455 0.09250 0.2740 1.8864 —1.0211 —5.0635 
7 0.01295 0.01874 0.03520 0.08491 0.2964 2.5820 —4.3869 
&mk 
0 0.04042 0.01096 0.004656 0.002573 0.001704 0.001295 0.001107 
eee 0.006723 0.01406 0.009831 0.007009 0.005339 0.004258 0.003932 
2 —0.012190 0.000358 0.003178 0.003382 0.003163 0.003068 0.002761 
3 —0.008555 —0.002515 0.000083 0.001020 0.001329 0.001415 0.001431 
4 —0.005308 —0.002468 —0.000789 0.000027 0.000396 0.000561 0.000619 
5 —0.003139 —0.001878 —0.000852 —0.000283 0.000011 0.000158 0.000236 
6 —0.001806 —0.001203 —0.000635 —0.000288 —0.000097 0.000005 0.000054 
: § —0.000823 —0.000581 —0.000329 —0.000167 —0.000073 —0.000023 0.000002 
hk 
0 0.01799 0.003558 0.001187 0.000550 0.000321 0.000225 0.000183 
1 0.004672 0.008717 0.005328 0.003404 0.002380 0.001769 0.001593 
| 2 —0.008196 0.000265 0.002257 0.002269 0.002012 0.001863 0.001644 
3 —0.004972 —0.001817 0.000062 0.000749 0.000950 0.000987 0.000983 
j 4 —0.002715 —0.001686 —0.000583 0.000020 0.000294 0.000412 0.000450 
5 —0.001546 —0.001221 —0.000614 —0.000210 0.000008 0.000118 0.000175 
6 —0.000783 —0.000746 —0.000457 —0.000212 —0.000073 0.000004 0.000040 
} Yj —0.000343 —0.000356 —0.000217 —0.000123 —0.000055 —0.000018 0.000001 
| TABLE 2 oc” 
e 
n Y F G H 
0 0.3709 ee ee Sia 2°} 
| 0.1951 0.3505 0.07361 0.008829 0.003227 
| 0.3827 0.3169 —0.11257 0.007196 0.003154 o" . rT) 
0.5556 0.2792 —0.20497 0.005715 0.002687 : 
0.7071 0.2400 —0.29460 0.004657 0.002244 é i | | i | ie cae Gl 
0.8315 0.1998  —0.43062 0.003952 0.001912 ae a. ee Se ee OEY 
0.92389 0.1536  —0.71476 0.003504 0.001691 ae soe | | 
‘ 0.9808 0.09381 —1.5335 0.003267 0.001570 TS See ee SE Re.” as 
CONCLUSIONS 1. —| 
The additional angle of attack a* through which pe SS Sas ee. | ‘ee . 
each swept-back wing section has to be turned in order | 
(2-4. 











to preserve the original lift distribution, was found to 

depend only on that lift distribution and, of course, on Fic. 3. Additional angle of attack a* for angles of sweep back 
the angle of sweep back. Thus, if two straight wings A= 1 end 00". 

have the same lift distribution but otherwise differ in 


section characteristics, they will, if swept back, lead to REFERENCES 


the same a*. 
The additional angle of attack is positive near the U.S. Bureau of Air Commerce, Spanwise Air-Load Distribu- 
tion, Washington, 1938, ANC-1 (1). 


root of the wing and negative near the tip. This was : : / 
t . ; ? Crean, J. R., Monoplane Wings with Sweeb—Theoretical Cal- 
o be expected, since sweep back causes a decrease in : tage nctapiny. gp gage : ED 

‘ 7 : culation of the Spanwise Lift Distribution, Aircraft Engineering, 
lift near the root and an increase near the tip and, to yo_ 10, pp. 245-247, 1938. 

compensate for that, the angle of attack has to be = s Weinig, F., Schiebende und gepfeilte Tragfltigel, Luftfahrt- 
changed accordingly. forschung, Vol. 14, pp. 45-54, 1937. 














